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Abstract 



The understanding of the physics of relativistic jets from active galactic nuclei (AGN) is 
still incomplete. A way to understand the different features of the AGN jets is to study 
it broadband spectra. In general, within the limits of present observations, AGN jets are 
observed in radio-to-X-ray energy band and they exhibit various intrinsic features such as 
knots. Particularly, the blazar jets which are pointed towards the observer, are observed in 
radio-to-7-ray and their radio maps exhibit internal jet structures. Moreover, the high energy 
emission from blazars show rapid variability. 

In this thesis, models have been developed to study the radiation emission processes from the 
knots of AGN jets as well as for blazar jets. A continuous injection plasma model is devel- 
oped to study the X-ray emission from the knots of sources 1 136-135, 1 150+497, 1354+195 
and 3C 371. The knot dynamics is then studied within the framework of internal shock 
model. In such a scenario, knots are formed due to the collision of two successive matter 
blobs emitted sporadically from the central engine of AGN. Shocks, generated in such col- 
lisions, accelerate electrons to relativistic energies. These electrons subsequently emit radi- 
ation via synchrotron and/or inverse Compton processes in the radio-to-X-ray energy range. 
The study of M87 knots involves a two zone model where the electrons with a power-law 
distribution are further accelerated. The synchrotron emission from these energetic electrons 
is then used to explain the observed spectrum. 

Regarding the blazar jets, the limb-brightening feature observed in the radio maps of the BL 
Lac object MKN501 is studied considering shear acceleration of electrons at the boundary 
of the jet. This interpretation does not require a large viewing angle of the jet as demanded 
by the earlier models and is consistent with the constraints obtained from very high energy 
studies. In case of MKN421, the dependence of temporal behaviour of radiation emission 
on the particle acceleration mechanism has been studied within the framework of two zone 
model. 
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Introduction 



1.1 Active Galactic Nuclei (AGN) 

The nucleus of a galaxy with luminosity > 10^"^ ergs s^^ and exceeding the overall luminos- 
ity of the entire host galaxy is known as Active Galactic Nucleus. Galaxies which host an 
AGN are called as active galaxies. 

1.1.1 Historical Background 

AGN was first observed in optical by Fath HI in 1908 at Lick Observatory. His aim was to 
test the claim of that time that the spectra of the spiral nebulae are continuous, consistent 
with a collection of stars. The continuous spectrum and absorption lines he observed for 
most of the sources suggested the presence of an unresolved collection of stars. However, 
for the nebula NGC 1068, the spectrum was composite showing emission and absorption 
lines. Later a higher quality and better resolution spectrum of NGC 1068 obtained by Slipher 
||2| at Lowell Observatory confirmed the presence of emission and absorption lines. Slipher 
also observed the emission lines are broad and spread over a substantial range of frequencies. 
Seyfert jSl was the first to do a systematic study of such galaxies with strong nuclear emission 
and broad lines. He attributed the broadness of emission lines to Doppler shifts and obtained 
a maximum velocity spread of 8500 km s^^ for the hydrogen lines of NGC 3516 and NGC 
7469. This class of objects with broad emission lines and a luminous core was later called 
as Seyfert galaxies (the most numerous type of AGN known). 

Though the observation of AGN began as early as the beginning of the twentieth century. 
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a major study of these objects started only after the advent of radio astronomy initiated by 
Karl Jansky [4], a radio engineer working at Bell Telephone Laboratories. Later Grote Reber 
flU, a radio engineer working on radio astronomy during his spare time, published a map 
of the radio sky at 160 MHz using his 31 feet reflector in his backyard. The map included 
the bright source located in the constellation Cygnus which we know at present as an AGN 
and commonly called as Cygnus A. In 1951, Smith [|6l used radio interferometry and ob- 
tained accurate positions of the radio sources Taurus A, Virgo A, Cygnus A and Cassiopeia 
A, which led Bade and Minkowski to identify the latter two in the optical band. The 
observed emission lines of these sources were very similar to that of Seyfert galaxies and 
radio sources with this characteristic feature in their spectrum are now identified as radio 
galaxies. There was also a considerable progress in the study of the radio source structures 
during this period. Jennison and Das Gupta [i8| used radio interferometry to study the struc- 
ture of Cygnus A which showed two bright components separated by ~ L5 arcmin. This 
morphology was observed to be common for extragalactic radio sources. Subsequent obser- 
vations at GHz frequencies provided better angular resolution and revealed the presence of 
a compact nucleus often associated with extended radio sources. In addition to these radio 
galaxies, some star-like radio sources were also detected. Their spectra were continuous with 
broad emission lines and no absorption lines. These peculiar stars were then called as "radio 
stars". Many attempts were made to understand them as an old novae or white dwarfs ([|9| 
and references therein) until Schmidt [UOll identified the lines of 3C 273 as nebular emission 
lines with redshift z = 0.158 and Greenstein and Matthews [[TTll identified emission lines in 
3C 48 with redshift z = 0.367 (the largest redshift known at that time). It was clear then 
that radio stars are highly luminous sources situated at large distances. These radio stars are 
today called as "quasi-stellar radio sources" or simply "quasars" (a term coined by Chinese- 
born U.S. astrophysicist Hong-Yee Chiu [fT2|). Seyfert galaxies, radio galaxies, quasars and 
objects of similar type are all collectively termed as Active Galactic Nuclei (AGN). Quasars 
are the distant AGN for which the host galaxy cannot be resolved. 

Janksy [fT3l suggested that the origin of the observed radio emission may be (a) from the 
stars or (b) secondary emission from the atmosphere due to the interaction of high energy 
particles emitted by the stars or (c) thermal emission from the interstellar dust. Former two 
options suggested maximum radio intensity from the direction of the sun since it being the 
closest star. However such an excess from the direction of the sun was not observed and 
hence these options are not viable. Meanwhile, Whipple and Greenstein [fT4l calculated the 
interstellar dust temperatures and it was found that they are too low to produce the radio 



Chapter 1 . Introduction 



3 



intensity observed by Jansky. Reber [fTSl suggested the free-free emission by ionised gas in 
the interstellar medium as a plausible mechanism for the observed radio emission. However, 
Henyey and Keenan [fT6l and Townes [jTTl showed that even this cannot reconcile the re- 
quired Jansky's brightness temperature. These studies thereby ruled out the possibility of the 
radio emission from the interstellar dust. Finally in 1950, Kiepenheuer [fT8l explained the 
galactic radio background in terms of synchrotron radiation by cosmic rays in the galactic 
magnetic field. By the end of 1950's, the synchrotron theory got accepted for explaining 
the radio emission from the extragalactic radio sources. Later the power-law' nature of the 
spectra and high degree of polarisation observed also supported this theory. 

In 1959, Burbidge Hl9l used synchrotron theory and found that the minimum energy content 
of the radio sources is extremely large ~ 10^° ergs. Burbidge attributed this large energy to 
a chain of supernovae explosions occurring in a tight packed set of stars at the nuclear region 
of the galaxy. Such a situation is plausible since the nuclear region is expected to be much 
denser than the average density of the galaxy. In 1962, Hoyle and Fowler [|20l suggested that 
such close packed stars can as well form a single super massive object. Accordingly they 
proposed the existence of a super massive star at the nuclear region of the galaxy and conjec- 
tured that it would evolve into a 'super- supernovae' providing the required energy. However 
the thermonuclear process is less efficient in converting mass into energy compared with the 
gravitational one for the masses of this order. Again, Hoyle and Fowler [.21] proposed in 
their pioneering paper in Nature in early 1963, that the energy of the radio sources was of 
gravitational origin derived from the slow contraction of a super massive object under its 
own strong gravitational field. This theory of a collapsed super massive object powering the 
AGN is widely accepted at present. We refer to this object as Central Engine hereafter. 

Today AGN are detected up to 7-ray energies and their structures are studied even at X-ray 
energies (see Chapter [3]). AGN are first detected in 7-rays by satellite based experiments 
followed by ground based experiments using atmospheric Cherenkov techniques. 3C 273 
was the first AGN detected in 7-ray (MeV) by Cos-B ^ in 1978 Later EGRET ^ op- 
erating at MeV-GeV 7-ray energy range detected around 60 AGN during 1991-2000 J23l. 
Around this period the ground based atmospheric Cherenkov telescopes (ACT) were opera- 

'We define a power-law flux as cx v^"^ (ergs cm^^ s^^ Hz^^) where v is the observed photon frequency 
and a the power-law spectral index. 

^ Cos-B was an European Space Research Orgarrisation satellite mission to study 7-ray sources launched by 
NASA. 

"'Energetic Gamma-Ray Experiment Telescope (EGRET) is a satellite borne 7-ray telescope onboard Comp- 
ton Gamma-Ray Observatory (CGRO) launched by NASA. 
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tional with the detection of MKN 421 in TeV 7-ray by Whipple ^ in 1992 [|24ll . Later imaging 
atmospheric Cherenkov telescopes (lACT) (e.g. HEGRA ^ TACTIC ^ CAT \ CANGAROO 
^ etc.) and wavefront sampling Cherenkov telescopes (e.g. CELESTE ^, PACT STACEE 
etc.) detected 6 AGN at TeV 7-ray energies [|25l . All AGN detected in 7-ray energies 
during this period are blazars (a class of AGN with a jet of matter flowing at relativistic 
speed towards the observer ( §1.1.4 and ^1.1.3)) except the radio galaxy Cen A which was 
detected by EGRET. Presently the satellite based experiment Fermi^^ operating at GeV-TeV 
7-ray energies detected around 600 AGN till now which include blazars and radio galaxies 
[|26ll . The second generation lACT namely MAGIC^^, HESS and VERITAS operating 
at GeV-TeV 7-ray energies detected around 40 AGN which again include blazars and radio 
galaxies'^. The number of AGN detected in 7-ray energy range are likely to increase with 
the help of these experiments and the future experiments (MACE'^ and CTA^^ ). Many of the 
AGN detected by Fermi and second generation lACT along with simultaneous observation 
at radio, optical and X-ray energies cannot be explained with our present understanding and 
the problems regarding AGN physics are still open [27J. 



"^Whipple is located at the Fred Lawrence Whipple Observatory in Southern Arizona, USA. 
^High-Energy-Gamma-Ray Astronomy (HEGRA) telescope is located at Roque de los Muchachos Obser- 
vatory on La Palma. 

^TeV Atomospheric Cherenkov Telescope with Imaging Camera (TACTIC) is located at Mt. Abu, Ra- 
jasthan, INDIA. 

^Cherenkov Array at Themis (CAT) is located at Themis, France 

^Collaboration of Australia and Nippon (Japan) for a GAmma Ray Observatory in the Outback (CANGA- 
ROO) is located at Woomera, Australia. 

'^CErenkov Low Energy Sampling and Timing Experiment (CELESTE) is located at Themis, France 
'"Pachmarhi Array of Cherenkov Telescopes (PACT) is located at Pachmarhi, Madhya Pradesh, INDIA 
"Solar Tower Atmospheric Cherenkov Effect Experiment (STACEE) is located near Albuquerque, New 
Mexico. 

^^Fermi Gamma-ray Space Telescope is a satellite based 7-ray telescope. The mission is a joint venture of 
NASA, the United States Department of Energy, and government agencies in France, Germany, Italy, Japan, 
and Sweden. 

'^Major Atmospheric Gamma-ray Imaging Cherenkov Telescope (MAGIC) is a system of two lACT situated 
at the Roque de los Muchachos Observatory, La Palma 

''^High Energy Stereoscopic System (HESS) is a stereoscopic lACT located at Namibia. 

'^Very Energetic Radiation Imaging Telescope Array System (VERITAS) is an array of four lACT located at 
Fred Lawrence Whipple Observatory in southern Arizona, USA. 

'^http://www.mppmu. mpg.de/ rwagner/sources/ 

'^Major Atmospheric Cerenkov Telescope Experiment (MACE) is an upcoming stereoscopic lACT at Hanle, 
India. 

'^Cherenkov Telescope Array (CTA) is a proposed open observatory and will consist of two arrays of lACT 
with one array at northern hemisphere and a second array at southern hemisphere. 
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Fig. 1.1: Very Large Array (VLA) map of the AGN 3C 175 at 4.9 GHz. M is the core, C and O are 
the hotspots and D to L are the knots. Figure reproduced from Bridle et al. ||28]| 

1.1.2 Morphology 

AGN comes in a variety of morphological structures and sizes. At one extreme we find an 
unresolved compact luminous core and in the other end we have, complex structures extend- 
ing up to hundreds of kiloparsec (kpc, Ipc = 3.0857 x 10^*^ cm). However one can develop a 
primary morphological feature of an AGN based on the commonly observed characteristics 
(Figure [TT]). 

• Core: These are compact unresolved luminous radio components coinciding with the 
nucleus of the associated galaxy (if resolved). They mostly have a spectrum which is 
a fiat power-law with an index a < 0.5. Cores are found in almost all quasars and in 
nearly 80 per cent of all radio galaxies. 

• Lobes: These are two extended regions of radio emission located on opposite sides of 
the galaxy or the nuclei. Lobes can be separated by several hundred kpc or in some 
extreme cases even up to a few megaparsec (Mpc) (e.g. 3C 236 has an overall size of 
4 Mpc). They have a steep power-law spectrum with index a > 0.5. Within the lobes 
one often finds local intensity maxima commonly called as hotspots. 



Chapter 1 . Introduction 



6 



AGN 




FRI 



FSRQ 



FSRQ 



SSRQ 



Fig. 1.2: AGN Classification 



Jets: These are narrow features that connect the compact core to the outer regions 
commonly referred to radio jets. They extend from pc scale to kpc scales. Radio jets 
often have bright regions along its length which are called as radio knots or simply 



knots. We shall consider the jets in detail in n.2 



1.1.3 Classification 



AGN are broadly classified into two groups, namely radio loud and radio quiet based on the 
ratio of their radio luminosity at 5 GHz to the optical luminosity at the B-band. Convention- 
ally the demarcating value of this ratio is taken around L^gbz/ Lb ~ 10 and roughly 15-20% 
of AGN are radio loud. 



Radio Quiet AGN 

The commonly observed radio quiet AGN are the Seyfert galaxies which have a morphology 
of spiral galaxies. They are classified based on their optical/UV properties. 
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• Seyfert Type 1: These are the Seyfert galaxies having a bright star-like nucleus and 
emit a strong continuum from far infrared to the X-ray band. Their spectrum contain 
broad emission lines with a width of the order of few thousand kilometers per second. 

• Seyfert Type 2: These are the Seyfert galaxies with weak continuum and narrow emis- 
sion lines with a line width of the order of few hundred kilometers per second. 

Seyfert galaxies with properties intermediate between type 1 and type 2 are classified as 
types 1.2, 1.5 and so on. Although Seyfert galaxies are categorized as radio quiet sources 
there also exists radio loud Seyfert galaxies detected even up to 7-ray energies ll29l . 

Radio Loud AGN 

Radio loud AGN are classified based on their morphology and optical/UV properties. 

• Fanarojf-Riley I (FR I) radio galaxies: These are extended sources often having mor- 
phology with core, jet and lobes. Their jets are often symmetric and become fainter as 
one approaches the outer extreme. Due to this feature they are called as edge-darkened 
sources. Their spectrum contain narrow emission lines features. 

• Fanarojf-Riley II (FR II) radio galaxies: These are more luminous than FR I radio 
galaxies with knotty jet and lobes with hot spots. In contrast to FR I type, their intensity 
falls off towards the nucleus and hence they are called as edge-brightened sources. 
Jets of these types of sources are well coUimated and often one sided. Their spectrum 
contain narrow emission lines features. 

FR I and FR II galaxies are often elliptical and together are referred to as Narrow Line 
Radio Galaxies (NLRG). 

• Broad Line Radio Galaxies (BLRG): These sources have a continuum and emission 
lines resembling those from Seyfert 1 galaxies. 

• Radio Quasars: These sources are distinguished from BLRG based on the luminosity 
and have a luminous nucleus which outshines the light from the host galaxy. These 



sources are often observed to have one sided jets with superluminal knots (see U .2. 1 ). 
They are further classified into Steep Spectrum Radio Quasars (SSRQ) and Flat Spec- 
trum Radio Quasars (FSRQ) depending on their radio spectral index is either a > 0.5 
or a < 0.5. 
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• BL Lacs: These sources have strong nuclear continuum with high polarisation and 
rapid flux variability. Their continuum extends from radio to 7-ray energies with few 
of them detected at TeV energies [|25l . Their emission lines are absent or weak. 

BL Lacs and FSRQ are collectively called as Blazars since they share many common prop- 
erties like, rapid variability, high and variable polarisation, superluminal motion etc. 



1.1.4 Unification 

The idea of unifying the different classes of AGN as a variant of objects belonging to a 
single population emanated from the fact that they share many common observational fea- 
tures. The most successful theory in AGN unification is based on the Orientation hypothesis, 
which assumes that the observed differences between the different classes of AGN are due 
to different orientations of similar objects with respect to the line of sight [|30l[3T1l as shown 



in Figure 1.3 



The unification of Seyfert 1 and Seyfert 2 is based on the fact that broad emission lines are 
present in the polarised spectra of the latter. Hence the unified picture derived for these 
sources assume that the broad line emitting regions are located closer to the central engine 
and a "dusty torus" is wrapped around the central region. Consequently the broad emission 
fines are hidden by the torus when the source is viewed edge-on (i.e. when the angle between 
the line of sight and the axis of the torus is large). On the other hand, narrow line emitting 
regions are farther from the central engine and hence are not obscured. Also, the unified 
picture assumes the presence of hot electrons far from the central region which are believed 
to reflect the broad lines towards the observer when the source is viewed edge-on. These 
reflected lines will therefore be visible in the polarised spectrum. Considering the above 
mentioned picture of the source, the edge-on view will resemble as Seyfert 2 while the face- 
on view as Seyfert 1 [|30ll . Also the unified picture of AGN contains twin relativistic jets 
emanating from the central engine, normal to the obscuring torus. High luminosity sources. 
Quasars and FR II are considered to be similar objects with the jets of the former aligned 
close to the line of sight (angle < 14°) and thereby causing relativistic beaming. Similarly, 
BL Lacs are the aligned jet version of FR I (low luminosity sources) [[3T1l . 




1.2 Jets of Active Galactic Nuclei 

The bridge of radiation connecting the central compact source with the extended lobes of the 
AGN is called as a "Jet". They are the common feature observed in most of the AGN and 
are either two-sided or one-sided. They may extend from pc to kpc scales and are visible 
in radio, optical and X-rays. Jets are interpreted as conduits for the transport of energetic 
particles from the nucleus to the extended radio structures at relativistic velocities. 



1.2.1 Evidence for Relativistic flow 
Superluminal Motion 

For many AGN, jet components are observed to move with velocities greater than the ve- 
locity of light and this phenomenon is referred to superluminal motion [32|. Rees [133]| 
interpreted this phenomenon as a result of bulk relativistic flow at an angle close to the line 
of sight. To understand this, let us consider a blob ejected from a stationary source with 
a velocity v at an angle il) with respect to the line of sight of the observer as shown in the 
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Fig. 1.4: A radiating blob moving at a relativistic speed v at an angle ip with respect to the line of 
sight of the observer 



Figure |1.4[ Let the blob emit signals for a duration dte, measured from the frame of the 
stationary source. The distance travelled by the blob during this time is v dt^. Due to the 
inclined motion, the signal emitted by the blob at the end of the duration dt^, travel lesser 
distance compared to the one emitted at the beginning of dt^ for a distant observer. Hence 
the observer will measure this interval as 

dto = dte{^l cosip j (1.1) 

where c is the velocity of light. Since the projected distance travelled by the blob in the sky 
plane is v dt^. sin^, the apparent velocity measured by the observer will be 

V dte sin if) 

Va = 



dto 



cos-?/' 



(1.2) 



or 



X = (1.3) 

1 — P COS ip 



where = Va/c and (3 = v/c are dimensionless velocities. In Figure 1.5 we plot the 
dependence of with the viewing angle ip for different relativistic motion of the blob. It 
is evident from the plot that the apparent speed can exceed c for the blobs with relativistic 
velocities and moving closer to the line of sight. 
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Gamma-ray transparency 

High energy 7-rays can interact with low energy photons to produce electron-positron pairs. 
The cross-section for this process is maximum when 

e^e^r^2{mecy (1.4) 

where and are the energy of the high energy and the low energy photons respectively 
and rUe is the electron mass. For example, a 100 MeV photon can pair produce with 5 keV 
X-ray photon. In general, the transparency of a medium to a radiation is described by its 
optical depth. Optical depth of a medium corresponding to a radiation of frequency p is 
defined as 

s 

r,{s) = J a,{s') ds' (1.5) 

so 

where sq ^ si^ the path of light travel and a,, is the absorption coefficient. The medium is 

said to be optically thick (or opaque) when r > 1 and optically thin (or transparent) when 
T < 1. The optical depth r to pair production for a 7-ray of energy in a homogeneous 
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region of size R can be approximated as [|34l 




7; - 



Air Rc 



(1.6) 



where L{2ex) is the luminosity of the target photon at energy 26^ and ctt is the Thompson 
cross-section. If the power-law spectra of AGN and quasars in the X-ray region extend to 
much higher energies, then the source may be opaque to 7-rays if r(e^) > 1. This situation 
can be conveniently stated in terms of compactness parameter (/) defined as [f35l 



where L is the total luminosity produced in the region of size R. The opacity condition, 
r > 1, can now be translated in terms of compactness parameter as / > 60. Estimation 
of compactness parameter requires the knowledge of X-ray luminosity and the size of the 
emission region. The latter can be estimated from the observed time variability of the flux as 



where t^ar is the variability time-scale. The inferred values of / for blazars from the observed 
X-ray luminosities are much larger than 60. For 3C 279, / ~ 5000 and for PKS 0528-1-134, 
/ ~ 15000. However EGRET mission detected many blazars at energies greater than 100 
MeV. In order to observe the 7-rays of these energies, the actual luminosity L of the source 
in its proper frame must be much smaller than the observed value and the actual size should 
be larger than the inferred value in order to fulfill the condition / < 60. Indeed, this can 
happen if the source is moving at relativistic speed towards the observer. The relativistic 
beaming effects will enhance the luminosity as [36] 



R^c 



t 



var 



(1.8) 



(5^ Li 



'int 



(1.9) 



and the size of the emission region will be reduced as 

T-, C ^var 
n 

6 



(1.10) 
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where Lobs and L^t are the observed and the intrinsic luminosity and the relativistic Doppler 
factor 6 is given by 

^=Fn — I M (1-11^ 

r(i — p COS'?/') 

Here, F = ( 1 — ) ^ is the bulk Lorentz factor, (3 = v/cis the dimensionless velocity and 
^jJ is the angle between the line of sight of the observer and the direction of motion of the 
source. The compactness parameter / in such a case will reduce to 

l = S-'h^^^ (1.12) 

Dondi and Ghisellini [|37l estimated the minimum value of the Doppler factor, corresponding 
to the optical depth r = 1, for the 7-ray bright blazars detected by EGRET. Their inferred 
minimum values of 6 were spread within 1.3 to 11.3. 



1.2.2 AGN Jet Features 

Lobes 

Lobes are the regions where the jets terminate and release their energy and momentum into 
the ambient intergalactic medium. The radio spectrum of the lobes are steep with index 
a > 0.5 and the emission is usually polarised. This suggests the lobes to be an optically 
thin synchrotron sources driven by efficient cooling of relativistic particles. They contain 
enhanced emission regions known as hot spots which are often coUinear with the central 
core (see Figure [TT[ ). In low resolution radio maps of the lobes, the presence of hot spots 
will make it appear as an edge-brightened source. The radio spectrum of the hot spots are 
flatter than the lobe spectrum with the index in the range a ~ 0.5 — 1. This suggests the 
hot spots as the location where the jet hit the ambient medium and the bulk kinetic energy of 
the beam is converted into the random energy through the shock formed at the collision (see 



^2.1.2). The shock accelerated energetic particles diffuse from hotspot to the lobes thereby 



providing a continuous supply of energy. 
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Knots 

These are the bright regions observed along the length of the AGN jet (Figure \L6^. The 



emission from these regions are strongly polarised and the polarisation angle is often per- 
pendicular to the jet axis. Knots are seen in radio, optical and X-ray images of AGN jets and 
many show superluminal motion. The non-thermal (power-law) nature of the knot spectrum 
along with the observed polarisation suggests that the low-energy(radio-to-optical) emission 
must be synchrotron radiation emitted by relativistic distribution of charged particles cooling 
in a magnetic field. However, the X-ray emission from the knots can be an extension of syn- 
chrotron spectrum itself or the inverse Compton radiation emitted via scattering of soft target 
photons by relativistic electrons (see chapter[3]). The target photons for the inverse Compton 
scattering can be either synchrotron photons itself or the photons external to the jet. The 
enhanced brightness of the knots is often interpreted as a result of efficient acceleration of 
charged particles probably by a shock present in the jet. 



Jet Asymmetry 

One of the intriguing characteristics of the AGN jets is that they are often observed to be 
one sided (Figure [Tj] ). The straightforward way to interpret this asymmetry is to relate it 
to relativistic beaming, since there are enough evidences that the jets are relativistic. If we 
assume that the central engine produces two similar jets ejected in opposite directions, then 
the brightness of the jet moving towards the observer will be enhanced due to relativistic 
beaming. Whereas the counter-jet will be dimmed (or invisible) since it moves away from 
the observer. If Fadv is the flux of the advancing jet and F^ec is the flux of the receding one, 
then the jet/counter-jet flux ratio J can be written as [[32ll 



j^F^^fl±J3co^ (1.13) 



Free V 1 ~ COS ll) 

where (3 is the bulk velocity of the jet in units of c and i}) is the angle between the jet direction 
and the line of sight of the observer. If a is the power-law spectral index of the intrinsic jet 
flux, then the index p can be either 2 + « in case of a continuous jet flow or 3 + a for a 



moving isotropic source. Using equation ( 1.3 1, J can be expressed in terms of superluminal 
velocity {(5 a) as 



(1.14) 
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Fig. 1.6: Image of 3C 273 with knots. Radio image by MERLIN (left), optical image by Hubble Space 
Telescope (middle) and X-ray image by Chandra with optical contour overlaid. Figure 
reproduced from Marshall et al. |[38]| 
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Fig. 1.7: VLA image of 3C 308 at 4.9 GHz. Figure reproduced from Bridle et al. lEl 



Here 5 is the relativistic Doppler factor given by equation ( |1.11[ ). The fact that all quasars 
are detected with one sided jet suggests that their jets are directed close to the line of sight 
of the observer. For example, the predicted value of J for the quasar 3C 273 is 9 x 10^ [f32l . 

The interpretation of the one sided jet as an outcome of relativistic beaming effects can be 
tested by measuring the difference in polarisation between the lobe in jet and counter-jet 
side. The lobe situated at the counter-jet side should be less polarised than the lobe at the 
jet side since the emission from the farther lobe face more depolarising medium along the 
line of sight. Indeed, such a difference was observed by Garrington and Conway [|39l for 49 
sources from a sample of 69, supporting the above interpretation. 



Bent Jets 

Sky maps of most AGN jets show curvature and in particular this feature is prominent in the 
radio trails or head-tail^^ sources. In high resolution radio maps of these sources, one can 
see jets emerging from the nucleus and bent through a large angle POl . For example, the sky 
map of the radio loud quasar PKS 2136-1-141 show a jet which is bent by an angle ~ 210° 



and this is the largest bending angle seen till today PTII (Figure |1.8| ). The bending of the 
radio jet can be interpreted in many ways. A distortion in the jet shape can occur as a result 
of ram pressure associated with the motion of the host galaxy through a dense intracluster 
medium. Bent jets with reflection symmetry can be explained if the host galaxy is associated 



'''These are sources with a radio morphology consisting of a bright nucleus and a single faint jet streaming 
on one side. 



Chapter 1 . Introduction 



17 




6 -5 -10 -15 -20 -25 -30 -35 

Relative R.A. (marcsec) 



Fig. 1.8: Very Large Baseline Array (VLBA) images of PKS 2136+141. Figure reproduced from 
Savolainen et al. BTI 



with a companion. The jet curvature, in such a case is the result of the acceleration of the 
parent galaxy introduced by its companion p2| . Jets of the radio source can also appear as 
an inversion symmetry if the jet precesses. In this case though the jet matter follows a linear 
path, it will appear curved in the sky plane due to precession. Alternatively, the jet can also 
be deflected in spite of the kinetic motion of the host galaxy or by geometrical effects as 
explained above. As the jet moves through the surrounding medium, it gets disrupted when 
the medium is an extremely dense intracluster gas. Pressure stratification of the intracluster 
gas can also cause the jet to curve P3ll44ll . 



Limb-brightened Jets 

Another feature seen in some AGN jets is limb-brightening at pc scales. High resolution 
radio maps of these jets show the edges are brighter than the central spine of the jet (Figure 



1.9). Such features are commonly observed in few FR I sources and blazars [|45ll . This 
feature is usually explained by the "spine- sheath" model where the velocity at the jet spine 
is larger compared to the velocity at the boundary. Such a radial stratification of velocity 
across the jet arises when jet moves through the ambient medium and the viscosity involved 
causes a shear at the boundary. The difference in the flow velocity of the jet between the 
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Fig. 1.9: VLA image of M87 with limb-brightened feature (circled). Figure reproduced from Owen 
et al. 61 

spine and the boundary can cause differential Doppler boosting and this in turn can produce 
a limb-brightened jet. Alternatively, velocity stratification at the jet boundary can accelerate 
the particles via shear acceleration P6l l47l (see chapter |4]). These particles can then emit 
radiation via synchrotron process giving rise to a limb-brightened feature. Particles at the jet 
boundary can also be accelerated by turbulent waves initiated by the instabilities [|48l . 



Chapter 2 



Particle Acceleration and Radiative 
Processes 



The observed very high energy 7-ray emission (up to GeV-TeV energies) from blazars sug- 
gests the presence of extremely relativistic particles in AGN jets. Also the power-law spec- 
tra observed over a broadband starting from radio-to-7-ray indicates the jet emission to be 
dominated by non-thermal processes. The promising mechanism by which particles can be 
accelerated to very high energies in AGN jets is "Fermi mechanism" suggested by Enrico 
Fermi [jSOl to explain the power-law nature of the cosmic ray spectrum. 



2.1 Fermi Acceleration 

In 1949, Fermi [|50ll proposed that particles can be accelerated to high energies when they 
are scattered by magnetic irregularities (magnetic mirrors) associated with moving clouds in 
the interstellar medium. Consider a cloud moving with a velocity V along the x-axis in the 



observer's frame (Figure 2.1). Let us assume the cloud to be massive in comparison with the 
scattered particle and hence the centre of mass frame (CM-frame) is the frame of the cloud 
itself. The energy and the momentum of the particle before scattering in the CM-frame can 



be obtained using Lorentz transformation as [|5T 



E' = T,{E-Vp,) (2.1) 
p', = TJp,-^^ (2.2) 
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Fig. 2.1: Illustration of Fermi acceleration mechanism: Collision between a particle with velocity v 
and a massive cloud moving with velocity V. (a) Head-on collision; (b) Follow-on collision. 
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where E and E' are the energy of the particle in the observer's frame and CM-frame respec- 
tively and px and p'^ are the corresponding momenta along x-axis (we represent the quantities 
in CM-frame with prime). The Lorentz factor of the moving cloud Tc is given by 



-1/2 



(2.3) 



Since the CM-frame is the frame of scatterer itself, the energy and the x-component of mo- 
mentum after collision will be 



Ei = E' 



Px,s 



-Px 



(2.4) 
(2.5) 



Transforming the scattered energy back into the observer frame we get 



E, = T,{E'-Vp', 



E 1 + 



-2Vp, 



(2.6) 



where we have used equations (2.1 1 and (2.2). If we express the momentum along x-axis in 
terms of energy 



E 



(2.7) 



The scattered energy can then be written as 



Es = TtE 1 + 



2 Vvx 



(2.8) 



Here is the velocity of the particle along x-axis. For V <^ c, the change in the particle 



energy due to collision can be obtained using equations (2.8 1 as 



AE = E,-E 



E 



V 



2 Vvx 



(2.9) 



where we have retained the terms only up to second order in V/ c. If the particle arrives at an 
angle 9 with respect to the velocity of the scatterer as shown in Figure 2.1[ then = —v cos 6 
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and we get 



AE V 

— = 2 - [vcose + V] (2.10) 
E 



From equation (2.10) it is evident that the particle gains energy in head-on collisions. On 
the other hand, for follow-on collisions cos 9 is negative, and the particle loses energy. Since 
in the scatterer's frame the velocity of the approaching particle is more than the receding 
one due to addition of velocities, the probabihty for head-on collision is more than that of 
follow-on collision. Hence there will be a net gain in particle energy. 



2.1.1 Stochastic Acceleration - Second Order Fermi Mechanism 



Let us consider a situation where the scattering centers move randomly, then the net gain 



is obtained by averaging equation (2.10) over the angle 9. The acceleration process is then 
referred to stochastic acceleration. An example of stochastic acceleration process is tur- 
bulent acceleration, where particles are energized via scattering by moving magnetic inho- 
mogeneities associated with turbulence in a flowing fluid. Since the collision probability 
between the particle and the scatterer is anisotropic, we need to know the probability of a 
collision happening at an angle 9 to perform the averaging. If we assume f ~ c, then the rate 
of collision in the scatterer's frame will be greater by a factor rc[l + (l^/c) cos 9] [|52l . Hence 
the probabihty of collision at angle 9 will be proportional to this factor and the average will 
be 



2V cos 9 



2 fV 



3 V c 



+1 



/ /i [1 + (\//c)/i] rf/i 

/[I + (v/c)/i] d^i 

-1 



(2.11) 



where = cos6'. From equation (2.10), the average energy gain will then be 



/ AE 



8 /V^ 
3 \~c 



(2.12) 



Since the average increase in energy is of the order {V/c)'^, this process is called as second 
order Fermi acceleration mechanism. 
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2.1.2 Shock Acceleration - First Order Fermi Mechanism 

If we consider only head-on collisions in Fermi acceleration, then the dominant term in 



equation (2.10) will be of the order of (V/c) (for v ^ c) and the acceleration process is 
called as first order Fermi mechanism. Acceleration of charged particles at a shock front 
is an example of this mechanism. A shock is a discontinuity initiated by perturbation in a 
supersonic fluid flow. Fluid on either side of the shock will be in different state of equilibrium 
and are connected by the conservation equations. The strength of a shock is determined by 
its Mach number M defined as 

M = — (2.13) 

Cs 

where vi is the velocity of the upstream fluid with respect to the shock front and Cs is the local 
speed of sound in the upstream region. Here, we denote upstream as the fluid ahead of the 
shock front and downstream as the one behind. Strong shocks are the one with M ^ 1. Par- 
ticles in a magneto hydrodynamic fluid (plasma) are scattered by magnetic inhomogeneities 
associated with turbulence in the flowing fluid. In presence of a shock in the fluid they get 
energised by crossing the shock front from upstream to downstream or vice versa. 

Let us consider the case of non-relativistic strong shock where the shock velocity U ^ Cg 
and U <^ c. In the shock frame the plasma will pass through it with an upstream velocity 
f i(= U) and the downstream velocity V2. From the mass conservation we get 

Pivi = P2V2 (2.14) 

where pi and p2 are the mass density of the upstream and downstream plasma. For a strong 
shock in the limit of M — )• 00 we can write [|53l 



Pi = Till (2.15) 

where 7^ is the ratio of the specific heats at constant pressure and volume. For a fully ionised 
gas we have 7^ = | and hence ^ = 4 and V2 = \U . Hence the plasma on either side of 
the shock (upstream or downstream) will see the plasma approaching from the other side of 
the shock (downstream or upstream) with velocity |f/. In the proper frame of the plasma 
(upstream or downstream), the particle distribution is isotropic due to scattering. Let us 
consider a particle of energy E with x-component momentum in the upstream plasma. 
We have chosen a coordinate system where x axis normal to the shock front. The energy of 
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the particle E' in the frame of downstream plasma will be 

E' = E+p^V (2.16) 
where V = jU is the velocity of the downstream plasma with respect to the upstream 



4 

E 
c 



plasma. Since the particles are relativistic, we can write p ^ E/c and ^ ^ cos 6. Here 6 
is the angle between the particle momentum and the shock normal. Hence the particle enters 
the downstream with an energy increment 

AE = E cos9 (2.17) 

The probability of the particle crossing the shock front and entering into the downstream 
within the angle interval 9 to 9 + d9 is proportional to sin9 cos 9 d9. Hence the average 
increase in energy in crossing the shock once is 



AE\ 2 /V 



E 3 V c 



(2.18) 



A similar situation happens for the particle crossing the shock front from downstream to 
upstream plasma. Hence the average fractional energy gain in making one round trip is 



4 [V 



E 3 V c 



(2.19) 



Thus a particle in the vicinity of a shock is scattered by magnetic inhomogeneities and gets 
accelerated to higher energy by crossing the shock front multiple number of times. Also, 
since the average energy gain is of the order of iy j c), shock acceleration is efficient than the 
stochastic acceleration. 



2.1.3 Acceleration at Shear Layers 

A velocity shear can arise in a plasma flow if it passes through a viscous medium. Parti- 
cles can then get accelerated when they are scattered between different velocity layers by 
magnetic irregularities [|46l (Figure |2.2| ). Let us consider the case of a gradual shear, where 
the mean free path of the particle is much smaller than the transverse width of the sheared 
velocity layer. Suppose the flow is non-relativistic and we choose a reference frame in which 
the local fluid is at rest and the flow velocity is along z-axis, U = [/^(x)^^. If r is the 
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mean scattering time, then the distance travelled by the particle along x-axis before getting 
scattered will be 

Sx= — coseT (2.20) 
m 

where pi is the momentum of the particle, m is its mass and 9 is the angle between the 
particle momentum and x-axis. The change in the fluid velocity due to this displacement 
will be 

Su=(^^^6xe, (2.21) 

where (^) is the shear velocity gradient. The momentum of the particle p2 relative to the 
local fluid will then be 

pi = pi \ I + 2 sm 6^ cos H ^ (2.22) 

V Pi Pi J 
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where is the angle between y-axis and the projection of the particle momentum in yz -plane. 

2 

Since E = the average fractional energy gain due to shear acceleration will be 



AE 



oc 



dx 



(2.23) 



The first order term in equation (2.22 ) will get cancelled because for every forward scattering 
of particle between two velocity layers there will be an equivalent reverse scattering with 
negative Su. 



2.1.4 Accelerated particle distribution 

The Fermi acceleration process leads to a power-law distribution of particles [[55l . Consider 
after each encounter the particle energy increases by a factor ^. i.e AE = ^E. Let us assume 
a particle with initial energy Eg enters the acceleration region. After n encounters the particle 
energy will be 

E = Eo{l + 0" (2.24) 

Also, let the probability of escape from the acceleration region after an encounter be Pgsc- 
Hence, after n encounters the probability of the particle to remain in the acceleration region 
will be (1 — Pesc)"- Thus the number of particles in the acceleration region with energies 
greater than E will be 



N{> E) oc ^(1 -P, 



esc) 



'I - p Y 

' (2.25) 



P 



esc 



From equation (2.24), n can be written as 



HE/E,) 

n = , , ^^ (2.26) 
ln(l + 



Substituting equation (2.26) in equation (2.25 ), we get 
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where the index p is given by 



in(r-f^) 



2.2 Radiation Emission Mechanisms 
2.2.1 Blackbody Radiation 

The radiation emitted by a distribution of charged particle which is in thermal equilibrium 
with itself as well as the emitted radiation is called as blackbody radiation. The specific 
energy density of a blackbody spectrum is given by Planck spectrum 

M^) = 1 ergs cm-=^ sr-1 (2.29) 

where u is the observed photon frequency, n the solid angle, h the Planck constant, fc^ the 
Boltzmann constant and T is the blackbody temperature. The peak photon frequency at 
which the energy density is maximum for a given temperature T is given by Wien's law. 



J^peak = 2.82{kB/h)T. Integrating equation (2.29) over the entire photon frequency and the 



solid angle we get the Stefan-Boltzmann law for the energy density of a blackbody spectrum 

u{T)=aT^ (2.30) 
where a = Aass/c and ass is the Stefan-Boltzmann constant. 



2.2.2 Bremsstrahlung Radiation 

When a charged particle moves in a Coulomb field it gets accelerated and emits electromag- 
netic spectrum. This radiation is known as Bremsstrahlung radiation. Let us consider the 
motion of an electron in the Coulomb field of an ion. We will assume small angle scattering 
where the deviation of the electron path from a straight line is negligible. The variation in the 
dipole moment of the electron-ion system gives rise to dipole radiation and it can be shown 
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that the total energy emitted per frequency by the electron is [|5T1l 

1 0, h<^v/uj 



(2.31) 



where u = 2'ku is the angular frequency of the emitted photon, h is the impact parameter, 
Ze is the charge of ion, v the electron velocity and nie and e are the electron mass and its 
charge respectively. The total power emitted per frequency per volume for a medium with 
ion density rij and electron density rie assuming a fixed electron speed v will be [ISTI 

dW levre^ . , , ,^ 

— TTT-r = — 1= neUiZ gff{v,u) (2.32) 

dudVdt sVSc^mlv ^^^^ ' ' 

where g// is known as Gaunt factor which is a function of energy of the electron and the fre- 
quency of the emitted photon. For a thermal distribution, the particles follow a Maxwellian 
velocity distribution and we obtain the power per frequency per volume as rtSTI 

dW 2^71 / 27r 



dudVdt 3 rup \3km, 



where g// is the velocity averaged Gaunt factor. From equation (2.33) we find that the 



bremsstrahlung emission from a thermal particle distribution gives rise to a flat spectrum 



with an exponential cutoff at about hu ~ kT (Figure 2.3 1. 



2.2.3 Synchrotron Radiation 

The radiation emitted due to the helical motion of a relativistic charged particle along a 
magnetic field is know as synchrotron radiation [l5n[56ll . The radiative power emitted by a 
relativistic electron with velocity /3(= v/c) moving at an angle a (pitch angle) with respect 
to an uniform magnetic field B is given by [|5T1l 

Psyn = lrlc^^ (3^ sin^ a (2.34) 

where re = e^/nie (? is the classical electron radius and 7 is the Lorentz factor of the elec- 
tron, 7 = (1 — For an isotropic distribution of mono energetic electrons we need to 
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Fig. 2.3: Intensity spectrum (arbitrary units) of bremsstrahlung radiation due to a thermal electron 
distribution at temperature 100 keV. 



average equation ( |2.34| ) over all angles and we get 



B 



(2.35) 



where ax is the Thomson cross section and Ub = B'^/8n is the magnetic field energy den- 
sity. Equation ( |2.35| ) is also equal to the average energy lost by an electron via synchrotron 
process. 

The synchrotron spectrum emitted by an electron of energy 7 nie (? moving with pitch angle 
a can be written as [|57l 



^3 e^B sin 



01 „ y 



nip 



(2.36) 



where v is the frequency of the emitted photon and 

36^72 



47r me c 



sma 



(2.37) 
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Fig. 2.4: Synchrotron power function 



is the critical frequency. The synchrotron power function F{x) is defined as 



oo 



(2.38) 



where K^/^ is the modified Bessel function of order 5/3. The shape of the spectrum is 



determined by F{x) with a peak located at ^ 0.29(z//z/c) (Figure 2.4). Alternatively one can 



write the single particle emission spectrum using equation (2.35) as [|5& 

4 

where is a function of v and 7 satisfying the relation 

M^f) du = 1 



(2.39) 



(2.40) 



Considering the shape of the spectrum (Figure 



2.4 1, one can approximate the function as a 
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Dirac delta function (to the orders of unity) 

0,(7)^5(z/-7V) (2.41) 

where z/^ = e5/27rme c is the Larmor frequency. 

For an isotropic power-law electron distribution given by 

N{'y) = k'J~P ^min<l<lmax (2A2) 

the radiation energy emitted per second per frequency for 'jmax ^ 7mm can be shown as lISTll 

a{p)iy-^P-^>/^ (2.43) 



dudt mf,c^ \47rmeC 

where a{p) is a function of particle spectral index. Hence the emitted synchrotron spectrum 
is a power-law with index {p — l)/2. 

Moreover, a charged particle in a magnetic field besides emitting synchrotron radiation, can 
absorb a photon and get energized. This absorption process is called as Synchrotron Self 
Absorption. Also a photon can induce a charged particle to emit in a direction and at a 
frequency of the photon itself {stimulated emission or negative absorption). The absorption 
coefficient for synchrotron self absorption process is given by [|58l 



1 



nrngu'^J 7 (^2 _ ]^^i/2 



,2 _ 1)1/2 ^ 'W ^syn{y,t^}\ 



(2.44) 



The specific intensity (1^) of the synchrotron radiation can then be found using the radiative 
transfer equation as [|56l 

I, = S,{1 - e-^") (2.45) 

where we have assumed a source with uniform properties and no background illumination. 
Here S,y(= iu/>^u) is the synchrotron source function and is the optical depth defined over 
a distance s as 

T^= I ds' (2.46) 
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The synchrotron emissivity for an isotropic emission is given by 



1 

477 



(2.47) 



For a power-law distribution of particle (equation (2.42)), the synchrotron self absorption 
coefficient will be [[56ll 



^ ^(p+2)/2 ^-(p+4)/2 



(2.48) 



Hence it may be possible for a source to be optically thick (r^, > 1) at low frequencies but 



optically thin (r^, < 1) at high frequencies. From equations (2.39 1, (2.47 1 and (2.48 ) one can 



find that the source function S,, for the optically thick region will be a power-law of the form 



5/2 



(2.49) 



The total synchrotron spectrum will then be a broken power-law with spectrum changing 
from u^/'^ at lower frequencies to //^(p^i)/^ at high frequencies. The frequency at which the 
index changes is called as synchrotron self absorption frequency. For AGN jet emission this 
frequency is observed to be within a range of few gigahertz. 



2.2.4 Inverse Compton Radiation 

Scattering of low energy electrons by high energy (or hard) photons is called as Compton 
scattering and the reverse process where high energy electrons scatter off low energy (or 
soft) photons is called as inverse Compton scattering. The spectrum obtained due to the 
scattering of soft photons by relativistic electrons is called as inverse Compton spectrum. If 
the energy of the incident photon in electron's rest frame is much smaller than the electron 
rest mass energy, then one can ignore the recoil of electron and the scattering process leave 
the photon energy unchanged in electron's rest frame. In such case the scattering process is 
described by Thomson cross section with differential cross section given by 

^^.^(l + cos^*) ,2.50) 

Here is the classical electron radius and 6 is the angle between the incident and the scat- 
tered photon directions. On the other hand if the recoil of the electron becomes considerable. 
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then the scattering cross section is described by Klein-Nishina cross section. The differential 
cross section in this case wiU include the quantum effects and is given by 



dar ri ez f t e 



e 



sm^e] (2.51) 



where e and are the energies of the incident and the scattered photon. These conditions can 
be expressed in terms of e and the Lorentz factor of the relativistic electron 7 as, 7e ^ nie 
for scattering in Thomson regime and 76 ^ mg for scattering in Klein-Nishina regime 
||5ni571 . Also, in the Thomson limit the scattered photon energy (e^) can be shown as ||571 

es^l^e (2.52) 

Hence for large 7 the photon energy gain is very large. Nevertheless this gain is quite small 
compared with the electron energy and the electron loses only a small fraction of energy in 
each scattering. However in Klein-Nishina regime the electron loses almost its entire energy 
to the photon in a single scattering. Hence the scattered photon energy in this case will be 

e, ^7mec2 (2.53) 

The radiative power emitted (or the power lost by an electron) due to inverse Compton scat- 
tering of an isotropic soft photon distribution can be shown as [1511 

Pcom = \l3^1^caTUph (2.54) 
where Uph is the energy density of the soft target photon distribution and the scattering is 



assumed to be in Thomson regime. Comparison of equation (2.54) with the power emitted 



due to synchrotron emission (equation (2.35)) one finds 



^ = ^ (2.55) 

-'com ^ ph 

When the scattering happens in Klein-Nishina regime the power lost by an electron can be 
computed using 



com,kn 




Chapter 2. Particle Acceleration and Radiative Processes 



34 



where is the emission rate of the scattered photon per frequency given by [|57l 



dN 2tt Tg c n(e) de 



dtdes 7 

X 



(2.57) 



Here n{e) is the number density of the soft target photons and the quantities and q are 
defined as 



4e7 



T^{-fmec'^ - es) 



(2.58) 
(2.59) 



For a power- law distribution of particle given by equation ( 2.42 1 the emitted photon spectrum 
when the scattering happens in Thomson regime will be [|57ll 



dW 

dtde^ 



+ Ap+ 11 



{p + 3)^p + l){p + 5) 



,-(p-l)/2 



X / e^P-^'^/^ n{e) de 



(2.60) 



From equations (2.43 1 and ( 2.60| ) we find that the spectrum emitted by both synchrotron and 
inverse Compton process can be represented by a power- law with same index (p — l)/2. In 
case of extreme Klein-Nishina limit the emitted spectrum will be 



dW 
dtde. 



2^|P+l g-p 



X 



nr^ck [rrie c ) 



- n{e) fin ^ + C{p) 



(2.61) 



where C{p) is a parameter of order unity [|57l . Thus we see the inverse Compton spectrum 



in extreme Klein-Nishina limit is much steeper than that of Thomson limit (equation (2.60 1) 



Synchrotron Self Compton 

In many astrophysical systems, the high energy emission is explained as a result of the in- 
verse Compton scattering of synchrotron photons. Here, the same electron population which 
is responsible for the synchrotron emission will scatter off these photons to higher energies. 
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This process is commonly referred as synchrotron self-Compton (or SSC) mechanism. For 
example, the ratio of the X-ray to TeV 7-ray fluxes obtained during simultaneous observa- 
tion of BL Lac objects, both in quiescent and flaring state, can be explained in the context of 
SSC process [|59ll . 

Let us assume that the radiation emitting plasma of AGN jet be confined in a spherical 
region with tangled magnetic field. This plasma moves down the jet at relativistic speed. 
The distribution of electrons in the emission region is assumed to be a power-law described 
by equation ( 2.42[ ). Therefore, the optically thin synchrotron spectrum will be a power-law 



with index a = (p — l)/2 (equation (2.43)). The SSC flux at the photon energy can then 



be predicted from the observed synchrotron flux as [|60l [Ml 

F^^^(e,) ^ d{a) ^^-2(2"+3) ^-(3a+5) ^-a ^1^^ '^"""'^ (2.62) 

where 9fi is the angular size of the source, is the synchrotron self absorption frequency, 
F^^" is the synchrotron flux at frequency z/^, Umax is synchrotron high frequency cutoff 
corresponding to the high energy cutoff in the particle spectrum, z is the redshift of the 
source and S is the Doppler factor of the jet. Here d{a) is a function depending only on a 
and has values rf(0.25) = 130, d{0.50) = 43, rf(0.75) = 18 and rf(l.OO) = 9.1. If we assume 
the X-ray emission from AGN jets as due to SSC process, then by comparing the predicted 
SSC flux at X-ray energy with the observed flux at that energy one can estimate the Doppler 
factor of the jet [|32ll . 



External Compton 

In external Compton mechanism, the photons which are produced outside the emission re- 
gion are scattered off to high energies by inverse Compton process. A general case assuming 
the scattering of an isotropic distribution of soft target photons to high energies by a distri- 



bution of relativistic electrons can be studied using equations (2.60) and/or (2.61 ). However 
in case of AGN jet these equations are invalid. Here as the emission region moves down 
the jet at relativistic speed, the plasma see an anisotropic distribution of target photons due 
to Doppler boosting. Consider the case of a spherical emission region moving down the 
jet with bulk Lorentz factor T. Let the inclination angle of the jet to the line of sight of 
the observer be Oq. For simplicity let us assume the external target photon distribution be 
monochromatic and isotropic in the frame of the central source. The energy density of these 
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target photons in the emission region frame will then be ~ V^u^^^, where m*,^ is the energy 
density of the isotropic external radiation field. The target photon distribution in the frame 
of the emission region will peak at energy e' ~ Fe*, where e* is the energy of the external 
photon distribution. The resultant inverse Compton emissivity jc{e, at photon energy e 
emitted in a direction ^l, due to a power-law distribution of particles described by equation 



( |Z42| ), is given by L62j| 

Ue, n) ^ ^^1%^ [ r (1 + ^) (4) " erg^ '^n^"' '''' (2.63) 

one \e / 

where a = {p — l)/2 and fi = cos 9z, with 9^ being the angle between the emitted photon 
and the jet axis. Since the scattering is assumed to happen in Thomson regime, equation 



(2.63 1 is valid only for the scattered photon energies satisfying the relation 



7^ • < < 7^ (2.64) 

imin — Tf*^! -|- y(x) — imax ■/ 

The quantity r(l + /i) can be written in terms of the viewing angle 9o and the Doppler factor 
of the jet 5 as 

r(l + /i) = 5^ — T- (2-65) 

1 + v/c 

where v is the velocity of the emission region along the jet. 



2.3 Equipartition Magnetic Field 



Consider a spherical source of volume V with a power-law distribution of electrons described 



by equation (2.42) cooling in a magnetic field B. The total energy of the electrons will be 



U, = Vm,^ / 7iV(7)rf7 



\ In 



2-p 



Imax ^min 



(2.66) 
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The total synchrotron luminosity of the source using equation ( |2.35| ) (assuming /3 ~ 1) will 
be 

'ymax 

L = V [ Psyn{l)N{-f)dj (2.67) 



ill-al - 7™D (2.68) 



syn \ 

AaTcUskV 



3(3-p) 

We can write 7mm and 7max in terms of characteristic synchrotron photon frequency z/„ 



and Vrnax using equation (241 1 as 7,,„i„ = {vmin/^LY''^ and 7max = {^'max/^LY''^- Then the 



ratio of total energy of the electrons to the synchrotron luminosity will be 

m A 



(2.69) 



where A is a constant that depends only on the particle spectral index p. If the source 
had other particles like hadrons along with electrons then the total particle energy will be 
Up = aUe, where a > 1. Then the total energy of the source will be 

Utot = Up + Ub (2.70) 
(^AL ^B^ ^^^^^ 

The magnetic field -Bmm for which the total energy of the system is minimum can be obtained 
by solving 



tot 



OB 



= (2.72) 



and we get 



Bmin = I -^^] (2-V3) 



V 



On the other hand, the equipartition magnetic field B^g obtained from the relation Up = Ub 
will be 
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Equations ( |2.73[ ) and ( |2.74[ ) differ by an factor less than 10 percent and the total energy 
corresponding to B^q and Bmin will be 



or 



UtotiBeq) 



Utot{Beq) ~ 



(r)2 
8? 

7 



and 



in 



(2.75) 
(2.76) 



7 V B^or, 



Utot{Bmin) ~ 1-01 ?7tot (-Bmin) 



(2.77) 



Hence it is customary to use equipartition value for the magnetic field while modelling the 
sources to ensure a minimum energy condition lfT9l[34ll . 



2.4 Hadronic Processes 



If protons are accelerated to high energies they can lose their energy through synchrotron 
emission and hadronic interactions. The main hadronic interactions by which an energetic 
proton can lose its energy are the following: 



Bethe-Heitler process: 

p + hu ^ p + e'^ + e~ (2.78) 



where hu is a photon. The photon threshold energy in the rest frame of proton for the 
Bethe-Heitler process is the sum of rest mass energies of the electron and positron (i.e. 
1.022 MeV). The cross section for this process in case of an ultra relativistic proton 
(/3 ~ 1) can be expressed in terms of the photon momentum k' when 2 < k' < A 
as (61 

, 271 . ,/A;'-2\V 1 23 , 37 , 61 A 

where afg is the fine structure constant, r^. is the classical electron radius, Z is the 
charge of the ion in units of e and r] = {k' — 2)/ {k' + 2). For k' > A the cross section 
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can be approximated as [|63l 



9 27 



6 In A;' - - + - In^ 2k' 
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27.512 



(2.80) 



proton-proton collision: 



p + p — )■ X + ^ TTj where X — hadrons 



(2.81) 



Here m is the multiplicity of secondary pions. The threshold energy for this reaction 
is Eth = 2m^c^(l + niT^/Amp) ^ 280 MeV, where and nip are the masses of 
the 7r°-meson and the proton. When the incident proton is in the GeV to TeV energy 
region, the total cross section can be approximated by [|64ll 



app{Ep) ^ 30 [0.95 + 0.06 In (Efci„/lGeV)] mbarn (2.82) 

where Ep is the initial proton energy and E^m = Ep — nipC^, and E^m > 1 GeV. Here 
it is assumed that cTpp = at lower energies. 

• photo-meson process: 

m 

p + hu ^ X + ^ vTj where X — > hadrons (2.83) 

1=1 

The cross section for the photo-meson process increases starting from the threshold 
energy of the photons Eth = 150 MeV (in the rest frame of protons) reaching their 
maximum value ~ 3 x 10"^*^ cm^ at ~ 300 - 400 MeV, and then decrease [|65l . 
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Decay modes for the pions produced in these reactions are as follows 

7r° ^ 27 (2.84) 
7r+(7r-) ^ + (2.85) 

/i+(/i-) ^ e^{e-) + Ue{iye) + JyM (2.86) 

The decay products, 7-rays and pairs, can then initiate an electromagnetic cascade by causing 
the production of further pairs and 7-rays. 



2.5 Emission Models 



The radio-to-UV/X-ray radiation from AGN jets are generally attributed to synchrotron 
emission due to cooling of relativistic non-thermal electrons in a magnetic field. However 
there are two different approaches concerning the high energy emission viz. leptonic model 
and hadronic model. High energy radiation ranging from MeV to TeV energies are gen- 
erally observed from blazars. Hence these models are discussed in the context of these 
sources. In leptonic models, the high energy radiation will be dominated by the inverse 
Compton emission from the same ultra relativistic electrons producing synchrotron radia- 
tion ll66l l67l [68l [69l . On the other hand, in hadronic models, the high energy radiation is 
mainly due to pair cascades initiated by the interaction of relativistic protons with photons 
and proton synchrotron radiation [rTOHTTll . 



2.5.1 Leptonic Models 

Leptonic models assume the electrons in the jets are accelerated to ultra relativistic velocities 
via Fermi acceleration process. Whereas protons are not sufficiently accelerated and their 
energies remain lower than the threshold energy required to initiate the hadronic interactions. 
Some models assume the hadrons are cold and mainly provide the inertia required for the jet 
to reach up to kpc/Mpc scales [|72l|73l. The accelerated electrons beside emitting synchrotron 
radiation also scatter off soft target photons to hard X-ray and 7-ray energies by inverse 
Compton process. The possible target photons for this process are the synchrotron photons 
produced within the jet (SSC) ( §2.2.4 ) and/or the external photons entering into the jet. The 



sources of external photons which can play an important role in explaining the high energy 
radiation in AGN jet are 
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• the radiation from the accretion disk around the central massive object 

• the reprocessed accretion disk radiation from the broad line emitting region 

• the infra-red photons from the dusty torus. 

Some parameters of the leptonic model can be constrained from the relativistic Doppler 
boosting required for the high energy radiation to be transparent against the pair production 
opacity with soft photons ( §1.2.1 1. However this effect may be non-negligible at very high 



energies and the resultant spectrum will be hard. Also, the synchrotron radiation from the 
secondary electrons may become important [[7411 . Moreover, the detection of subluminal 
velocities (/3app < 1) in the sub-pc scale jets of few TeV blazars suggest that the relativistic 
jets of these sources decelerate. The varying Doppler factor due to this deceleration will have 
a significant impact on the observed properties of the blazars [1751 . 

In simplistic approaches, the underlying electron distribution is either a single or broken 
power-law with index/indices inferred from the observed photon spectral index/indices. A 
reasonable estimate of the parameters can then be obtained from the observations. For ex- 
ample, the underlying magnetic field can be estimated considering equipartition between the 
electrons and the magnetic field energy densities ( §2.3 1. As discussed earlier, the Doppler 



factor can be estimated from the observed superluminal motion of the knots or using 7- 



ray transparency (^1.2.1). Also, from the measured variability time-scale, the size of the 
emission region can be constrained. While these simplistic models have been successful in 
reproducing the spectrum of AGN jets, they lack a self-consistent basis for the shape of the 
electron distribution. 

A more realistic approach consists of the solution of a kinetic equation involving acceleration 
of the particles and radiative as well as non-radiative cooling mechanisms [[761 . One zone 
models assume the observed emission as a result of efficient cooling of non-thermal electrons 
from a region with tangled magnetic field. The evolution of the particle distribution A^(7, t) 
in this region can be described in its simplest form by the kinetic equation as 

^^^^1^ + ^ [P(7, t)iV(7, t)] = g(7, t) (2.87) 

where 7 is the Lorentz factor of the electron, P(7, t) is the energy loss rate and Q{'y, t) is 
the injection rate of non-thermal particles. The injection can be a single burst of non-thermal 
particles injected at time t = {one-time injection models) or a continuous injection of non- 
thermal particles. If the losses are mainly due to synchrotron and inverse Compton processes. 
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then for a power-law distribution of particles similar to the one given in equation ( |2.42[ ), the 
resultant photon spectrum will be a power- law with index {p — l)/2. However since the 
energy loss rate is proportional to 7^ for these processes (^2.2.3 and ^2.2.4 ), the high energy 
particles cool more efficiently than the low energy ones. This leads to a depletion of high 
energy particles in case of one time injection models and gives rise to a time-dependent 
high energy cut off in the non-thermal particle distribution. The emitted photon spectrum 
will then exponentially decrease at high energies corresponding to this cut off energy in the 
particle spectrum. Alternatively one can infer the age of the emission region by translating 
this exponentially decreasing feature in the photon spectrum to the high energy cut off in 
the particle distribution. On the other hand, in the continuous injection models, the depleted 
high energy electrons are continuously replenished. This gives rise to a broken power-law 
particle distribution with a break at energy for which the cooling time-scale is equal to the 
age of the emission region. The resultant photon spectrum will then be a broken power-law 
instead of one with an exponential cut off. 

Two zone models are more involved than one zone models where acceleration of particles 
are also considered along with the cooling mechanisms. According to these models particles 
are accelerated to relativistic energies in an acceleration region. These high energy particles 
are then injected into a cooling region where they lose most of their energies by radiative and 
non-radiative processes. The evolution of the particles are governed by the kinetic equations 
corresponding to acceleration region and cooling region. These equations can be written in 
their simplest form as 

an(7, t) ^^^^^^^^ ^ ^^^^^ ^^^^ ^ ^ ^^^^ ^2.88) 



Ot O'J t 



esc 



where the equation (2.88 1 governs the evolution in the acceleration region and the equation 



(2.89) in the cooling region. Here 72(7, t) and A^(7, t) are the particle distribution in the ac- 
celeration region and the cooling region, PARil, t) and Pcr{i-, t) are the respective energy 
loss rates, 'jacc is the particle acceleration rate, tesc is the particle escape time-scale in accel- 
eration region and Q{'j,t) is the particle injection rate. The injection into the acceleration 
region can be mono energetic electrons or a residual particle distribution of an earlier accel- 
eration process. If the acceleration happens at a shock front then the acceleration rate jacc 
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can be approximated as 

iacc ~ 7^ (2.90) 

where t^cc is the acceleration time-scale and can be estimated from the theory of diffusive 
shock acceleration [77] . The index of the particle spectrum in acceleration region is governed 
by the acceleration and escape time- scales. Whereas the maximum energy to which the 
particles can be accelerated is determined by the acceleration and cooling time- scales. One 
zone and two zone models involve many parameters even in their simplistic form and are 
often cannot be constrained. However these models can throw light on the underlying physics 
of the source and its predictions compared with future observations will help us to understand 
these sources better. 

Sambruna et al. [|78l used a one-time injection model to explain the broadband emission 
from the knots of several AGN detected in radio, optical and X-ray. The spectra modelled by 
them predicts an exponential cut off at optical/UV or X-ray energies. Whereas Liu & Shen 
ITTQII proposed a two zone model to explain the X-ray emission from the knots of M87 (a 
nearby FRI radio galaxy). One zone models failed to reproduce the observed flux and/or the 
spectral index from the knots of this source. 

In case of blazars the picture is different. One-zone models are used by various authors to 
explain the emission from blazars llSOllSTI . These models assume a homogeneous distribu- 
tion of particles and magnetic field throughout the emission region. The broadband spectra 
of the blazars are successfully explained by this model. However the time lags observed 
between the flares at different frequencies cannot be explained under this model. Also this 
model requires the light travel time to be shorter than the synchrotron cooling time- scales 
in order to satisfy the homogeneity of the particle distribution. This introduces a constraint 
on the size of the emission region. Two zone models offers more insight into these sources 
ll82l[83l . Kirk et al. [[82ll solved the kinetic equation considering the temporal as well as the 
spatial variation of the particle distribution in the cooling region. They showed that the time 
lags between the flares at different frequencies can be an outcome of the difference between 
cooling and acceleration time- scales. 
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2.5.2 Hadronic Models 



In hadronic models, the protons are accelerated along with electrons to ultra-relativistic ener- 
gies and cool off mainly through proton synchrotron and photo-meson interactions (synchrot- 
ron-proton blazar model) [ItTII . The dominant channels for the photo-meson process are [.84,1 

p + 7 — )■ 7r° + p 
p + 7 — )■ vr^ + n 



The pions decay as shown in equations (2.84), ( 2.85[ ) and (2.86) and decay products initiate 
an electromagnetic cascade. Hadronic models explain the high energy emission as a result of 
these cascades and proton synchrotron emission. The type of the resulting cascade spectrum 
depends upon the compactness parameter / which measures the optical depth with respect 
to pair creation (equation ( |1.7[ )). For emission region with small compactness parameter 
/ ^ 60 the electromagnetic cascade terminates after few generations and hence the photon 
luminosity is concentrated at high energies. On the other hand, as / increases, more and more 
generations shift power towards lower energies. 

The target photons for photo-meson process can be the electron synchrotron radiation and/or 
the external photons. Miicke & Protheroe [85] using Monte Carlo technique simulated the 
proton interactions and the subsequent cascades. They considered the co-acceleration of 
protons along with electrons while the synchrotron emission from the latter is responsible 
for the low energy emission from blazars. These photons serve as target photons for the p7 
interactions. They showed that the cascades initiated by the 7r° decay and n"^ decay generate 
a featureless 7-ray spectra. In contrast, the proton synchrotron cascades and /i^ synchrotron 
cascades produce a two-component 7-ray spectrum commonly observed in flaring blazars. In 
general, direct proton and synchrotron radiation is mainly responsible for the high energy 
bump in blazars, whereas the low energy bump is dominated by synchrotron radiation from 
the primary electrons, with a contribution from the secondary electrons [1861 . 

Unlike leptonic models, hadronic blazar models result in neutrino emission through the pro- 
duction and decay of charged mesons (equations ( |2.85[ ) and (2.86)). Another important 
source of high energy neutrinos is the production and decay of charged kaons. In the case 
of p7 interactions positively charged kaons are produced [[87l[88l . They decay into muons 
and direct high energy muon-neutrinos. These muon-neutrinos will not have suffered energy 
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losses through tt"^ and fi"^ synchrotron radiation unlike the ones originating from n'^ and /i^ 
decay. Therefore they appear as an excess in comparison to the remaining neutrino flavors 
at the high energy end of the emerging neutrino spectrum. Detection of predicted neutrino 
spectrum play an important role in validating the hadronic models. 

Investigation of time-dependent hadronic models is very difficult because of time consum- 
ing Monte-Carlo cascade simulations. Also it is difficult to reconcile their rapid variability 
observed in blazars (< 1 hour) with the radiative cooling time-scales of protons ff89l . 

3C279 




V [Hz] 

Fig. 2.5: Spectral fits to the spectral energy distribution of 3C279 using a leptonic external-Compton 
model (solid (red)); leptonic SSC model (short-dashed (red)); hadronic model with electron 
synchyrotron photons as targets for the photo-meson process (dot-dashed (maroon)) and 
hadronic model with electron synchrotron -i- external photons as targets for photo-meson 
process (long-dashed (maroon)). Figure reproduced from Bottcher et al. 11271 

Bottcher et al. ll27l studied the simultaneous multi wavelength observation of the very high 
energy (VHE) blazar 3C279 using the leptonic and hadronic models. Leptonic one zone 
model requires unrealistic parameters to explain the observed spectrum. Whereas, the hadronic 
synchrotron-proton blazar model is able to fit the broadband spectrum successfully. They 
also considered the contribution from the external target photons for hadronic interactions 
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in order to reduce the energy loss time-scale. In Figure 2.5 we show their fit to the spectral 



energy distribution of 3C279 based on leptonic and hadronic models. 



2.6 Aim of the Thesis 

Despite the availability of enormous amount of information about AGN by virtue of high 
resolution and high sensitivity experiments at present, there exist a large amount of uncer- 
tainties regarding the physics of various observed features [|90l |9T1 |92l |93l. In the work 
presented in this thesis we shall attempt to understand certain features of AGN jets in view 
of the recent observations. 

• In Chapter |3| we describe our works to interpret the underlying physics of the AGN 
knots. The on-board X-ray satellite Chandra^ studied the knots of several AGN for 
which radio and optical informations are already available [|94l |95l 1781 |96ll . The X- 
ray emission from these knots can be either due to synchrotron emission or inverse 
Compton emission [|9l|95l|96l|971|78l|9l. If the X-ray flux lies below the extrapola- 
tion of radio-to-optical flux then the synchrotron origin is plausible else the emission 
may be due to inverse Compton process. Sambruna et al. [|78l used one zone model 
to explain the X-ray emission from these knots. This model gives rise to a spectrum 
with a time-dependent exponential high-frequency cutoff. Hence for the knots with 
synchrotron origin of X-ray, their model predicted an exponentially decreasing steep 
spectrum at this energy. However the photon spectral index measured from the short 
duration observations of the knots at X-ray energies contradicts this prediction [98]. 
On the other hand, the acceleration process may exist for a longer duration such that a 
continuous injection of non-thermal particles may be viable. We model the observed 
radio-optical-X-ray spectra of the knots of the AGN 1136-135, 1150-1-497, 1354-1-195 
and 3C 371 by using a continuous injection plasma model. We assume the knot to 
be a uniform expanding sphere with continuous injection of non-thermal particles. 
The electron distribution and resultant radiation spectrum is computed by taking into 
account synchrotron cooling, inverse Compton scattering of cosmic microwave back- 
ground and adiabatic cooling due to the expansion of the sphere. The continuous 
injection of particles will generate a break in the electron distribution at energy where 
the cooling time-scale is equal to the age of the emission region (cooling break). The 

' Chandra is a satellite borne X-ray telescope launched by NASA. 
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energy at which the index changes is time dependent, and synchrotron/adiabatic cool- 
ing is dominant for the particles with energy greater than this break energy. We show 
that this model can successfully reproduce the observed spectrum from the knots of 
these sources. 

We also interpret the knots as an outcome of an internal shock which we discuss in 
§3.2[ An internal shock interpretation for the knots was first suggested by Rees ||99ll 
to explain the knots of M87. We study the jet dynamics and the viability of the above 
mentioned continuous injection model using a simple internal shock model. The cen- 
tral engine of the AGN is assumed to emits blobs of matter sporadically at relativistic 
speeds. In the process, the fast moving blobs will collide with the previously ejected 
slow moving ones thereby forming a shock. We implement the model for the knots 
of AGN and compute the time-evolution of the non-thermal particles produced. Also 
we compare the results obtained with the broadband fluxes from knots of several AGN 
jets and their observed positions. The motivation here is to find quantitative values of 
the model parameters by demanding that the model can self consistently explain the 
observation. 

It is also noted that the observed X-ray flux from the knots in the jet of the nearby 
galaxy M87 cannot be explained by considering simple one zone models involving 
continuous injection or one-time injection of non-thermal particles HIOOII . Perlman & 
Wilson mOOII proposed a modified CI model where the volume within which particle 
acceleration occurs is energy dependent. Using this phenomenological model they ob- 
served that the particle acceleration takes place in a larger fraction of the jet volume in 
the inner jet than the outer jet. Also particle acceleration region occupy a smaller frac- 
tion of the jet volume at higher energies. Liu & Shen [[791 proposed a two-zone model 
with the acceleration region and cooling region spatially separated. The advection of 
particles from the acceleration region to the cooling region introduces a break in the 
particle spectrum which along with the cooling break in the cooling region produces a 
double broken power law. The synchrotron emission from such a particle distribution 
is used to fit the observed spectra. We propose an alternate two-zone model to explain 



the observed X-ray flux from these knots and which we discuss in ^: 3.3 Here we con- 
sider the acceleration of power-law distribution of electrons in an acceleration region 
which are then injected into a cooling region and cool via radiative processes. We show 
that one will obtain a broken power-law particle spectrum when the initial spectrum 
injected into the acceleration region is flatter than the characteristic spectrum of the 
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acceleration region. This particle distribution will then develop an additional cooling 
break depending upon the age of the knot while cooling in the cooling region. The 
synchrotron emission from the resultant particle distribution in the cooling region is 
used to reproduce the broadband spectrum of the knots of M87 jet. 

• In Chapter |4} we discuss a model to explain the observed limb-brightened structure 
of the blazar MKN 501 . The high resolution image of the nearby BL Lac object MKN 
501 in radio show a transverse jet structure with the edges being brighter than the 
central spine. Such a feature is commonly referred as "limb-brightened" structure. 
This feature is usually explained by the "spine-sheath" model where the velocity at the 
spine is larger than the velocity at the boundary. The limb-brightened structure can 
then be explained by considering the differential Doppler boosting for a proper combi- 
nation of the bulk Lorentz factors and the viewing angle. The viewing angle deduced 
for MKN 501 based on this model is > 15° HIOIH . However the high-energy studies 
of MKN 501 demand the viewing angle to be ~ 5°. Since the high-energy emission is 
originated from the inner part of the jet close to the nucleus, a possible bending of the 
jet was suggested by the earlier work HlOlll . We explain the observed limb-brightened 
structure of MKN 501 jet as an outcome of efficient particle acceleration process at 
the jet boundary. Here the particles can be accelerated by shear acceleration or turbu- 
lent acceleration. We deduce the required condition for the shear acceleration to be 
dominant over turbulent acceleration and discuss the diffusion of particles accelerated 
at the boundary into the jet medium. Also we derive the spectral index of the parti- 
cle distribution accelerated via shear acceleration process and turbulent acceleration 
process and show the observed index at the boundary of MKN 501 jet supports the 
former. The bending of the jet is not a requirement for this interpretation unlike the 
explanation based on differential Doppler boosting. 

• In Chapter |5| we study the temporal and spectral behaviour of the non-thermal emis- 
sion from blazars. The flux variability observed in blazars have been studied by several 
authors using one zone and two zone models [l80l [82l [T02l [83l . Kirk et al. [82J and 
Kusunose et al. [83] assumed a two zone model where particles are accelerated in a 
region presumably by a shock and escape into emission region where they lose their 
energy by radiative processes. Chiaberge & Ghisellini H 10211 explained the short time 
variability observed in MKN 421 by dividing the emission region into thin slices. We 
study the spectral and temporal behaviour of blazars by considering a two zone model 
under two different scenarios of acceleration process. Mono energetic particles are 
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assumed to be accelerated in the acceleration region and then injected into the cooling 
region where they lose most of their energy via radiative processes. The behaviour of 
the resultant spectrum is then studied for two cases of particle acceleration process in 
the acceleration region. In the first case, the rate of particle acceleration is assumed to 
be energy dependent and in the second case, it is independent of the same. The model 
is then applied on the BL Lac object MKN 421 and its flare characteristics are studied 
under these two cases. It is found that, detailed information about the temporal be- 
haviour of blazars can throw light on the underlying particle acceleration mechanism. 

• Finally in Chapter |6} we summarize the work presented in the thesis and discuss the 
possible future work. 



Chapter 3 



Emission Models and the Dynamics of 
AGN knots 



Knots in kpc scale jets of several AGN are recently studied in X-ray by the on-board X-ray 
satellite Chandra [|94l |95l 1781 |96l . Chandra due to its excellent spatial resolution is able 
to resolve bright X-ray knots and in most of the cases it coincides with their radio/optical 
counterparts 1195117811981 . The radio-to-optical emission from these knots are generally ac- 
cepted to be of synchrotron origin, whereas the X-ray emission could be due to synchrotron 
[1971 1781 l98l or inverse Compton processes depending on its radio-to-optical (aRo) and 
optical-to-X-ray index (aox) [|95l|78l|98l|96l|96ll93. If oro > aox, then the X-ray flux 
hes above the extrapolation of radio-to-optical flux and hence a single emission mechanism 
may not explain the observed fluxes. In such a case, the X-ray emission may be due to in- 
verse Compton process or it may arise from a different electron population other than the 
one responsible for the radio/optical emission. However, it should be noted here that such a 
spectrum can still be a resultant synchrotron radiation from an electron distribution modified 
due to inverse Compton scattering happening at extreme Klein- Nishina regime [I103II . The 
synchrotron spectrum in this case can satisfy the spectral indices requirement, a^o > olqx- 
On the other hand, if olro < ceox, the X-ray flux lies below the extrapolation of radio-to- 
optical flux and synchrotron origin of X-ray is acceptable [|96l I78l l98l . The synchrotron 
origin of X-rays for knots was strengthened in case of the knots of 3C271 since the alternate 
inverse Compton model would require exceptionally large Doppler factors [|98ll . 

When the X-ray emission can be attributed to the inverse Compton process, the possible 
choices of target photons are radio/optical synchrotron photons (SSC) Hi 0411 or external pho- 
tons. The source of external photons which can be dominant at kpc scale jet is the cosmic 
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microwave background' (IC/CMB) [|96l|95l|78l|9l. The SSC interpretation of X-ray emis- 
sion would require large jet powers and magnetic fields much lower than the equipartition 
values whereas IC/CMB requires relatively low jet power and near equipartition magnetic 
fields [1961. 



These possible radiative process identifications have to be associated with (and confirmed 
by) dynamical models regarding the origin and subsequent evolution of the radiating non- 
thermal particles. In many models, these non-thermal particles are assumed to be generated 
by a short duration acceleration process and the particle distribution is determined by radia- 
tive losses (one-time injection) [|95l |78l |96l |105[ l76l I106II . The high energy particles cool 
more efficiently ( |2.2.3 and {2.2.A) and hence get depleted in time. This give rise to a time- 



dependent high energy cut off in the non-thermal particle distribution. If the X-ray emission 
is attributed to synchrotron emission by these particles, then these models predict an expo- 
nentially decreasing X-ray spectrum [19511781196119811 . This can be translated to a high energy 
cutoff in the electron distribution which in turn gives an estimate of the age of the knot. 
These non-thermal electrons move with a bulk speed v c along the jet. Thus from the age 
of the knot, one can determine the location in the jet where the short duration acceleration 
process occurred. The distance of the knot from the central object and the short duration of 
acceleration (much less than the age of the knot) may naturally put strong constraints on any 
models of the acceleration process. Also the photon spectral slope measured during the short 
duration observations of 3C 371 is ax ~ 1.7 ± 0.4 which is in apparent contradiction to the 
predicted exponential X-ray spectrum by this model [[981 . Moreover, the model requires the 
coincidence that the age of the knot be equal to the time required for X-ray emitting elec- 
trons to cool. A larger survey of X-ray jets have to be sampled to confirm whether this is 
statistically plausible. 

On the other hand, it may also be possible that the acceleration process exists for a duration 
longer than the age of the knot, and hence, there is a continuous injection of non-thermal 
particles. The acceleration process may be due to internal shocks formed as a result of colli- 
sion between sporadically ejected relativistic blobs of matter from the central engine [[107[[ . 
In such a case the duration of acceleration can be roughly the cross-over time of the collided 
blobs. If we assume a typical blob size ~ kpc (which is roughly the size of the knots seen in 
radio) and its velocity ~ c, the acceleration duration can be as large as ~ lO^^s. In the next 
section we study the case where the knots of AGN are modelled as an spherical expanding 

'According to the Big Bang theory, cosmic microwave background (CMB) radiation is the relic radiation 
left over from the formation of the universe. 
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emission region with a continuous injection of non-thermal particles and discuss the results 



obtained. In ^: 3.2 we interpret the same considering an internal shock scenario and study its 
dynamic properties. Though one-time injection and continuous injection models can suc- 
cessfully reproduce the observed spectrum of knots of many AGN, they failed to explain the 
same for the knots of a nearby AGN, M87. The X-ray flux of M87 knots when compared 
with radio-to-optical flux suggests a synchrotron origin but its flux and/or spectral index can- 
not be reproduced by these simple models. In §3.3| we propose a modified synchrotron model 
to understand the broadband emission from the knots of M87. 



3.1 A Continuous Injection Plasma Model 

We consider the knot of the AGN as a plasma moving relativisticaUy along the jet with a 
bulk Lorentz factor T. In the rest frame, it is assumed that the plasma uniformly occupies an 
expanding sphere with radius R(t) = Ro + (3expCt, where Ro is the initial size of the sphere 
and (3expC is the expansion velocity. Initially at t = there are no non-thermal particles in 
the system. A continuous and constant particle injection rate for t > is assumed, with a 
power-law distribution of energy, 

Q{-f) d-f = Qo'y'P d-f for 7 > 7™^ (3.1) 

where 7 is the Lorentz factor of the electrons. The evolution of the total number of non- 
thermal particles in the system, A^(7, t), can be conveniently described by a kinetic equation 
of form 

+ ^ [^(7, t)iV(7, t)] = g(7) (3.2) 
Here P{'y,t) is the particle energy loss rate given by 

P(7, t) = -{is{t) + iic{t) + iA{t)) (3.3) 



where 7s (t), iicif) and 7a (i) are the cooling rates due to synchrotron, inverse Compton 
scattering of the cosmic microwave background (IC/CMB) radiation and adiabatic expansion 
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respectively. These cooling rates are given by 

, , 4 ctt B^(t) . 

75^=77 — (3.4) 
3 me c Sir 

licit) = ^ ^ asBTL,{z) 7^ (3.5) 

Mt) - (3.6) 

The former two cooling rates are associated with radiative losses and the latter is a non- 
radiative energy loss associated with the loss in internal energy used up for expansion [1761 . 
Here the evolving magnetic field is parameterized to be B{t) = Bo{R{t) / RoY"" and Tcmb{z) = 
2.73(1 + z) is the temperature of the CMB radiation at the redshift z of the source. Note that 
the time t and other quantities in the above equations are in the rest frame of the plasma. 



We solved the equation (3.2) numerically for N{'y,t) using the finite difference scheme 
described by Chang & Cooper HlOSlI and the resultant synchrotron and inverse Compton 
spectra are computed at an observing time t = to. Finally, the flux at the Earth is com- 
puted taking into account the Doppler boosting [f36ll . characterized by the Doppler factor 
5 = [r(l — (3 COS 9)]~^, where (3c is the bulk velocity of the plasma moving down the jet and 
9 is the angle between the jet and the line of sight of the observer^. 

While the total non-thermal particle distribution has to be computed numerically, a qualita- 
tive description is possible by comparing cooling time-scales with the observation time to. 
The cooling time-scale due to synchrotron and inverse Compton cooling at a given time t 
and Lorentz factor 7 is tc{t, 7) ~ 7/(75 + 7/c)- Then, 7^, defined as the 7 for which this 
cooling time-scale is equal to the observation time, tc{to, 7c) ~ to, becomes, 

^ me c t~^ 
^^"^[^ + 4nca5BTlJz)] 

The adiabatic cooling time-scale ta also turns out to be ^ to, since ta ~ R{t)/l3expC ~ to for 
Ro <^ R{to). Thus, the non-thermal particle distribution at time t = to can be divided into 
three distinct regions 11761 : 

1. In the regime 7 <C 7c, radiative cooling is not important and A^(7, to) ~ Qol'^io. The 



Here and everywhere else in this thesis Ho = 75 km s ^ Mpc ^ and go = 0.5 are adopted 
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corresponding spectral index for both synchrotron and inverse Compton emission are 

a = {p-l)/2. 

2. In the regime 7 3> 7c, either synchrotron or inverse Compton cooling is dominant 
and N{'y,to) oc 7~(p+^). The corresponding spectral index for both synchrotron and 
inverse Compton emission are a = p/2. 

3. In the regime 7 ~ 7c, either synchrotron or inverse Compton cooling as well as adia- 
batic cooling are important and the spectral slope is in the range (p — l)/2 > a > p/2. 

The computed spectrum depends on the following ten parameters: the observation time to, 
the magnetic field at the time of observation Bf = B{t = to), the magnetic field variation 
index m, the radius of the knot at the time of observation, Rf = R{t = to), the index p, 
the minimum Lorentz factor 7mm, the Doppler factor S, the bulk Lorentz factor T, expansion 
velocity (3expC and the normalization of the injection rate Qo- On the other hand, there are 
only three observational points, namely, the radio, optical and X-ray fluxes. Clearly, the 
parameters are under constrained, and it is not possible to extract meaningful quantitative 
estimates. However, the motivation here is to show that this model can explain the observed 
data with reasonable values of the above parameters. 



3.1.1 Results and Discussion 

We applied the above model to the knots of the AGN observed by Chandra along with the 



information available in radio and optical energies (Figure |3.1[ |3.2[ |3.3[ ). In Figure |3.4[ 
the computed spectra are compared with the data for different knots for four sources namely 
1136-135, 1150-1-497, 1354-1-195 and3C 371. The values of the parameters used are tabulated 



in Table 3.1 The injected power in non-thermal particles in the rest frame of the knot can be 



written as 



P.n, = I hm^c') Qil) dl = qo'^ 7jr'^ (3.8) 



p-2 

while the total jet power can be approximated as [17211 



Pj,t = vrP' /3ciUp + U, + Ub) (3.9) 
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where Up, Ue and U b are the energy densities of the protons, electrons and the magnetic field 
respectively. Here it has been assumed that the protons are cold and the number of protons is 
equal to the number of electrons. We found the jet power ranges from 10''® to 2 x 10"^® ergs 
s~^, while the injected power is generally three orders of magnitude lower. This means that 
the non-thermal acceleration process is inefficient and most of the jet power is expected to 
be carried to the lobes. The magnetic field Bf is nearly equal to the equipartition values. 

The X-ray emission for the knots in 3C371 and Knot A of 1 136-135 are identified as being 
due to synchrotron emission which is consistent with earlier works fTSl l98ll . However, in 
this case the predicted X-ray spectral index is ax = aji + 1/2 instead of being exponential. 
Note that this relation between the spectral indices is independent of the parameters used to 
fit the data. For the rest of the sources the X-ray emission is attributed to IC/CMB which is 
again consistent with the results obtained from the earlier works [|78l . However, for some of 
the sources the optical spectral index is now ao = an + 1/2 instead of being exponential. 

The present model can be confirmed (or ruled out) vis-a-vis one-time injection models, by 
future measurements of the radio a^, optical ao and X-ray ax spectral indices. In particular, 
the following cases are possible: 
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Relative R.A, (arcsec) Relative R,A. (arcsec) 



Fig. 3.2: Chandra image of 1150+497 and 1354+195 overlaid with VLA radio contours. Figure 
reproduced from Sambruna et al. fTSl 
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Fig. 3.3: Multi wavelength image of 3C371 in X-ray by Chandra{top), in optical by HST (middle) 
and in radio by Merlin (radio). Figure reproduced from Pesce et al. fl98l1 
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Fig. 3.4: The observed fluxes in radio, optical and X-ray compared with model spectrum using pa- 

I, while the 



rameters given in table 3.1 The data for 3C 371 is taken from Pesce et al. 
rest are taken from Sambruna et al. [78|. Triangles correspond to knot A and circles to 
knots. 
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Table 3.1: Parameters for model fittings 



v-.J V > Ul IVllWL 


Br 

(xlO^s G) 


limn 


F 


(xlO" s) 




p 


R 

Hexp 


P 

(ergs s"'-) 


p. , 

(ergs s"-*-) 


Bt IB 


1136-135A 


0.9 


2.0 


2.4 


0.2 


5 


5 


0.8 


44.7 


47.5 


0.4 


1136-135B 


4.0 


20.0 


2.9 


9 


5 


5 


0.1 


44.2 


47.9 


0.5 


11 50+497 A 


2.5 


30.0 


2.85 


9 


5 


3.5 


0.1 


44.2 


47.3 


0.4 


1150+497B 


4.3 


30.0 


3.3 


9 


5 


3.5 


0.1 


44.1 


47.3 


0.75 


1354+ 195 A 


1.7 


40.0 


3.0 


9 


3.5 


2 


0.1 


45.8 


48.4 


0.04 


1354+195B 


8.0 


25.0 


3.2 


9 


3.5 


2 


0.1 


44.7 


47.5 


0.63 


3C 371A 


1.3 


10.0 


2.4 


12 


3.5 


3.5 


0.1 


42.8 


46.4 


0.8 


3C 371B 


1.0 


10.0 


2.4 


1 


3.5 


3.5 


0.5 


43.5 


46.0 


0.9 



Columns: (1) Source and knot name taken from Pesce et al. |l98l for 3C371 and the rest from Sambruna et 
al. 1.78 J : (2) Magnetic field at the observation time, B ^ B{t = to)', (3) Minimum Lorentz factor jmin', (4) 
Power- law index of the injected non-thermal particle p; (5) Observation time to', (6) Doppler factor (5; (7) Bulk 
Lorentz factor F; (8) Velocity of expansion /S^xp (in units of c); (9) Log of the injected power Pinj', (10) Log of 
the total jet power Pjet', (10) Ratio of the magnetic field to the equipartition value. For all cases, the magnetic 
field variation index m and the size of the source at t = to is fixed at 1.5 and 5 x 10^^ cm, respectively. 

1. In the case aji ~ ax, the X-ray emission is probably due to IC/CMB. Both the 
continuous injection and one-time injection models are equally viable. 

2. In the case ax ~ ctij + 1/2, the X-ray emission would be due to synchrotron 
emission from electrons in the cooling dominated region. The continuous injection 
scenario will be favored in such case. 

3. In the case ax > «i? + 1/2, when the X-ray emission is exponentially decreasing, it 
should be attributed to the high energy cutoff in the electron distribution. The 
one-time injection scenario will be favored. 

4. In the case ax < olr, when the X-ray emission is exponentially increasing, it should 
be attributed to the low energy cutoff (7mm) in the electron distribution and the X-ray 
emission should be due to IC/CMB. Both the continuous injection and one-time 
injection models are equally viable. 

Similar arguments can be put forth for the optical spectral index ao as compared to the 
radio. It should be noted that for some older systems the one-time injection would be the 
natural scenario, while for younger systems the continuous injection would be more prob- 
able. The technique described above will be able to differentiate between the two, and a 
generic constraint on the acceleration time- scales and typical age of the knots may be ob- 
tained. A generic model where the injection rate decays in time may then be used to fit the 
observations. The measurement of spectral indices in different wave-lengths will also reduce 
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the number of unconstrained parameters in the model fitting, leading to reliable estimates of 
the system parameters. 

3.2 Internal Shock Interpretation 

The possibility of a continuous injection of non-thermal particles into the knots of AGN as 
well as its dynamical properties can be studied using a specific injection model. While there 
is no consensus on the origin of these non-thermal particles, one of the standard model is 
the internal shock scenario [|99l I107L Here the particles are energized by Fermi accelera- 
tion in shocks produced during the interaction of relativistically moving blobs ejected from 
the central engines with different speeds. A detailed description of the shock formation and 
subsequent electron acceleration is complicated and would require numerically difficult mag- 
neto hydrodynamic simulations. Moreover from the limited number of observables, which 
can be obtained from the featureless spectrum in two or three different energy bands, one 
may not be able to constrain the various assumptions and/or the initial conditions of such 
a detailed study. Nevertheless, a qualitative idea as to whether the internal shock model is 
consistent with the present observations (and if so, qualitative estimates of the model param- 
eters) would be desirable. Such an estimate would provide insight into the temporal behavior 
of the central engine. 

We implement an internal shock model with simplifying assumptions and compute the time 
evolution of the non-thermal particles produced. We also compare the results obtained with 
the broadband fluxes from knots of several AGN jets and their observed positions. The mo- 
tivation here is to find a consistent set of model parameters that can explain the observations 
and thereby make qualitative estimates of their values. Apart from the fluxes at different 
energy bands, the spectral indices in each band can also provide important diagnostic infor- 
mation about the nature of these sources. Hence, we have analyzed long (> 40 ks) Chandra 
observations of three AGN and present the constrains that were obtained on the X-ray spec- 
tral indices of the individual knots. 

3.2.1 Chandra X-ray Data Analysis 

Long-exposure Chandra observations of the sources PKS 1136-135, PKS 1150-1-497, and 
3C 371 were performed with the Advanced CCD Imaging Spectrometer (ACIS-S) with the 
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Table 3.2: Chandra Observations 



Source name Obs Id Exposure Knots i^o. 3-3.0 ctx 

(ks) (ergs cm~^ s~l) 



1136- 


135 


3973 


77.37 


A 
B 


0.63 
1.56 


1 r,^+1.51 

U.DO_Q 28 


1150 + 


497 


3974 


68.50 


A 
B 


3.15 
0.61 


0.66l°i« 

o.9ol?:g2 


3C371 




2959 


40.86 


A 
B 


3.32 
7.84 


1 43+0.85 

1 07+°-2^ 

l-^'-0.23 



Columns: (1) Source name; (2) Chandra Observation Id; (3) Exposure time; 
(4) Knots prominent in X-ray; (5) Flux in 0.3 — 3.0 keV energy band; (6) X-ray 
energy spectral index. 



source at the aim point of the S3 chip. The Observation ID (ObsID) and the exposure time 



of the observation are given in Table 3.2 Earlier shorter duration observations of these 
sources revealed two bright knots for each source, whose positions from the nucleus are 



given in Table 3.3 These longer duration observations allow for better constraint on the 
X-ray spectral indices of these knots. 

The data from Chandra X-ray observatory^ were analyzed using the Chandra data analysis 
software CIAO"^ and the latest calibration files were used to produce the spectrum. The X-ray 
counts from each individual knot was extracted using a circular region centered at the knot. 
The background was estimated from the counts obtained from same size regions located 
at the same distance from the nucleus but at different azimuth angles. The radius of the 
circular region was chosen to be 0.74" for the sources 1136 — 135 and 1150 + 497, while for 
3C371 a smaller radius of 0.6" was used. These sizes were chosen to minimize any possible 
contamination from the nucleus and/or the other knot. 

Spectral fits were undertaken on the data using the XSPEC package^ in C statistic mode, 
which is the appropriate statistic when the total counts are low. The flux and the energy 



spectral indices obtained are tabulated in Table 3.2 The quoted spectral indices are not very 
meaningful due to large errors. 

^http://cxc. harvard.edu/cda/ 
'*http://cxc.harv ard.edu/ciao/ 
^http://heasarc.gsfc. nasa.gov/docs/xanadu/xspec/ 
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Table 3.3: Observed Knot Features 



Source name 


Type 


z 


Knot 


Position 

(arcsec) 


OIRO 


aox 


Ref 


1136-135 


FSRQ 


0.554 


A 


4.5 


0.73 


0.83 


|78| 








B 


6.7 


1.04 


0.68 




1150 + 497 


FSRQ 


0.334 


A 


2.1 


0.99 


0.83 


m\ 








B 


4.3 


1.24 


0.44 




1354 + 195 


FSRQ 


0.720 


A 


1.7 


1.05 


0.6 


|78| 








B 


3.6 


1.15 


0.68 




3C273 


QSO 


0.158 


A 


13 


0.86 


0.61 


isi 








B 


15 


0.9 


0.73 




3C371 


Bl Lac 


0.051 


A 


1.7 


0.9 


1.28 


|98| 








B 


3.1 


0.76 


1.14 





Columns: (1) Source name; (2) Type of the source; (3) Redshift; (4) Knots 
prominent in X-ray (The nomenclature is same for all the knots as they are in 
the literature except for 3C371 where A and B are reversed); (5) Position of the 
knot; (6) Radio-to-Optical index; (7) Optical-to-X-ray index; (8) References: 
Sambruna et al. f78l, Sambruna et al. 1^, Pesce et al. f98l. 



3.2.2 Model 

In the internal shock model framework, temporal variations of the ejection process produces 
density fluctuations (moving with different velocities), which collide at some distance from 
the source to produce an observable knot. This distance will depend on the time-scale over 
which the variation takes place. In general the system will exhibit variations over a wide 
range of time-scales and knots like features would be produced at different distance scales. 
Here, we consider large-scale jets (with deprojected distances ^100 kpc) which are expected 
to arise from variability occurring on a corresponding large time-scale. Variations on smaller 
time-scales would produce knot structures on smaller distance scales, for example, pc scale 
or even smaller jets, which would be unresolved for these sources. These smaller time- 
scale variabilities will be smoothed out at large distances, and hence one expects the jet 
structure of these sources to be determined by variations over a single characteristic time- 
scale. To further simplify the model, we approximate the density and velocity fluctuations 
as two discrete blobs with equal masses. Mi = M2 = M, having Lorentz factors, Fi and r2 
that are ejected one after the other, from the central engine with a time delay of Ati2. The 
collision of the blobs is considered to be completely inelastic; i.e. the blobs coalesce and 
move as a single cloud, which is identified with the observed knot. From conservation of 
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momentum, the Lorentz factor of the knot is 



ri/3i + r2/32' ' 



1 (3.10) 



where f3i^2 = ^1,2/c are the velocities of the blobs. Since the collision is inelastic, a fraction 
of the bulk kinetic energy is dissipated. Denoting all quantities in the rest frame of the knot 
by subscript K, this dissipated energy AE^ can be estimated as 

AEK = [r,K + {r2K~2)]Mc' (3.11) 

The Lorentz factors of the blobs in the knot's rest frame T1K.2K = (1 ~ (^ik2k)^^^'^ 
computed using 

filK,2K = f'''~ (3.12) 
-L — Pi, 2 P 

where /3 = (1 — l/F^)^/^. The time-scale on which this energy will be dissipated can be 
approximated to be the crossing-over time of the two blobs, 

ToN,K ^ —5 (3.13) 

C (P2X - PlK) 

where, Axk is the average size of the two blobs in the rest frame of the knot. It is assumed 
that this dissipated bulk kinetic energy, AE^, gets converted efficiently to the energy of the 
non-thermal particles produced during the collision. The number of non-thermal particles 
injected per unit time into the knot is taken to be, 

QKij) d'y = A'y-P d-f for j > jmin, (3.14) 

where 7 is the Lorentz factor of the electrons, p is the particle index and A is the normaliza- 
tion constant given by. 



AEk (p-2) ^ (p_2) 
ToN,K rrieC^ 



A = ^ ^ ^ -f^^r') (3 15) 

m ^2 I mm 



Here the injection is assumed to be uniformly occurring for a time Ton,k- The cloud is 
assumed to be permeated with a tangled magnetic field, Bk- 



We used the kinetic equation (3.2) to study the evolution of the total number of non-thermal 
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particles in the system. However for the total loss rate we considered only the synchrotron 



and IC/CMB (equations (3.4) and (3.5)). Using equation (3.6), the adiabatic cooling time- 
scale can be written as ^ 7/7^1 = R{t) / {(3exp c). And for any given time of observation, 
to,K, the size of the blob will be R{to,K) = Ro + (^expCto,K- From the above one finds 
that ta will be always larger than to,K and hence the adiabatic cooling can be neglected. For 
to,K < ToN,K, the resultant particle distribution will be a broken power-law giving rise to a 



composite synchrotron spectrum with a spectral break (^3.1 ). 



The predicted spectrum and size of a knot depends on nine parameters, which are the mass 
of the blobs M, their average size Axk, the Lorentz factors Fi and the particle injection 
index p, the minimum Lorentz factor 7mm, the magnetic field Bq, the inclination angle of the 
jet 6, and the observation time to,K- 

From these parameters and the location of the knot in the sky plane, one can infer the time 
delay Ati2 between the ejection of the two blobs. The projected distance of the knot from 
the source S can be written as 

S = c{/3itc + l3to) sm9 (3.16) 

where to = ^to,K is the time of the observation after the formation of the knot in the source 
frame. The time elapsed tc, after the ejection of the first blob and the start of the collision, is 
given by, 

= "^^^"-^^^ (3.17) 
where Axi is the size of the first blob ^ Fi xAx^^/F. Thus Ati2 can be estimated using 



the above equation, where tc is given by equation (3.16), and it essentially depends on four 
parameters, 6, Fi, F2 and to^K- 

A total time, ttot, can be defined to be the time that has elapsed between the ejection of the 
first blob and the observation, ttot = tc + to- For two knots, A and B, the time difference 
between the ejection of their first blobs, t^^ is then 

= ttot - ttot - thT (3.18) 



where tiT is the light travel time difference between the two knots, which is approximated 
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to be 

cA cB 

tiT ~ , . (3.19) 

The power of the jet, can be defined in two different ways. The instantaneous power, which 
is the power when the system is active, can be defined to be the average energy of the blobs 
divided by the time-scale on which the blobs are ejected. This power can be estimated for 
each knot to be 

Mc^fTi + T2)/2 

Pins ~ — (3.20) 

On the other hand, the time averaged power of the jet can be defined as the typical energy 
ejected during active periods divided by the time-scale on which such activity occurs. For 
two knots, A and B, this can be approximated to be 

^ [MV(r^ + r^) + M^c^Tf + rf )]/2 

^ave ~ ^^^^ P.^i; 

3.2.3 Results and Discussion 

The model has been applied to those knots of kpc scale jets, which have been detected by 
Chandra and for which radio and optical data are available. This criterion was satisfied 



by the two brightest knots of the AGN: 1136 - 135 (Figure 3.1 ), 1150 + 497 (Figure 3.2), 



1354 + 195 (Figure 3^273 (Figure[r6]) and 3C371 (Figure[33]). In this work, the knot 
closer to the nucleus is referred to Knot A and the further one as Knot B. For these sources 
this nomenclature is same as in the literature [|78l[95l| except for 3C371, for which the farther 



one has been referred to Knot A [|98]| . For three of these sources, the X-ray spectral indices 



were constrained using long exposure observations as described in 3.2.1 The observed 



properties of the sources and the knots are tabulated in Table 3.2 and 3.3 



Figure 3.5 shows the observed radio, optical and X-ray fluxes of these knots along with the 
computed spectrum corresponding to model parameters that are given in Table |3.4[ Since 
the number of parameters is large as compared to the observables, a unique set of parameter 
values cannot be obtained. Two consistency checks have been imposed on the parameter 
values: that the number of non-thermal electrons which will be injected into the system, 
Nnth, is smaller than the total number of protons, A^^, and that the magnetic field, B should 
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10 15 20 25 

log i/(Hz) 

Fig. 3.5: The observed fluxes in radio, optical and X-ray compared with model spectrum using pa- 



rameters given in Table 3.4 Knot A fluxes are represented by filled triangles, and knot B 
fluxes are represented by filled circles. 
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Table 3.4: Internal Shock Model Parameters 



Source 


Knot 





Ti 


r2 




to,K 




^tot 


logM 


Imin 


P 


B 






(deg) 








(10"s) 


(IQi^s) 


(IQi^s) 


(g) 






(lO-^G) 


1136-135 


A 


11.5 


4.6 


5.4 


5.0 


0.25 


11.8 


11.6 


38.0 


2 


2.4 


1.1 




B 




4.1 


5.9 


5.0 


8.5 


13.1 


17.3 


36.4 


20 


2.9 


4.5 


1150+497 


A 


10.2 


2.5 


3.0 


2.8 


16.0 


1.2 


5.1 


37.5 


30 


2.8 


1.4 




B 




2.6 


3.5 


3.0 


6.2 


8.8 


10.0 


36.8 


30 


3.3 


4.0 


1354+195 


A 


8.21 


1.7 


2.3 


2.0 


9.0 


8.1 


8.1 


38.1 


37 


3.0 


1.7 




B 




1.7 


2.3 


2.0 


87 


15.0 


17.0 


37.1 


20 


3.2 


8.0 


3C 273 


A 


8.23 


2.2 


4.0 


3.1 


2.1 


27.0 


24.0 


36.8 


50 


2.7 


0.4 




B 




1.7 


2.3 


2.0 


8.0 


30.0 


32.0 


37.8 


80 


2.8 


0.6 


3C 371 


A 


15.8 


1.6 


2.4 


2.0 


5.9 


0.18 


1.4 


35.5 


20 


2.5 


1.0 




B 




2.0 


2.4 


2.2 


2.3 


0.39 


0.75 


36.8 


10 


2.4 


0.6 



Eight model parameters and derived quantities. The ninth parameter is /\xk = 5-0 x 10^^ cm for aU sources. 
Columns marked with an asterisk are derived quantities and not parameters. 

Columns: (1) Source name; (2) Knot; (3) Viewing angle; (4) Lorentz factor of the first blob; (5) Lorentz factor 
of the second blob; (6) Lorentz factor of the Knot; (7) Observation time; (8) Collision time; (9) Total time; (10) 
Mass of the blobs; (11) Minimum Lorentz factor of the particle injected into the knot; (12) Injected particle 
spectral index; (13) Magnetic field. 



Table 3.5: Knot/Jet Properties 



Source 


Knot 


Ati2 
(10"s) 


(10"s) 


B 

Be,. 


N„th 
Nk 


log Pins 
(ergs s-i) 


log Pave 
(ergs s-i) 


ToN,K 
(10"s) 


tO,K 
ToN,K 


D 

(kpc) 


1136-135 


A 


1.1 


3.1 


0.55 


0.88 


48.60 


48.17 


20.4 


0.01 


112.8 




B 


2.5 




0.74 


0.74 


46.65 




9.1 


0.93 


163.4 


1150+497 


A 


0.9 


5.0 


0.12 


0.13 


47.98 


47.31 


17.0 


0.94 


50.5 




B 


3.9 




0.65 


0.43 


46.63 




10.7 


0.58 


96.1 


1354+195 


A 


7.5 


18.0 


0.04 


0.38 


47.49 


47.15 


9.6 


0.94 


787 




B 


17.0 




0.64 


0.78 


46.10 




9.6 


0.90 


135.5 


3C 273 


A 


20.0 


34.0 


0.03 


0.80 


46.00 


46.63 


5.4 


0.39 


240.7 




B 


31.6 




0.02 


0.16 


46.60 




9.6 


0.84 


248.5 


3C 371 


A 


1.4 


1.5 


0.39 


0.88 


45.59 


46.98 


7.2 


0.83 


11.2 




B 


1.0 




0.28 


0.29 


47.14 




16.3 


0.14 


77 



Columns: (1) Source name; (2) Knot; (3) Time delay between the ejection of the blobs; (4) Time delay between 
the ejection of the first and the third blob; (5) Ratio of the magnetic field to equipartition magnetic field; (6) 
Ratio of non-thermal electrons to the total number of protons; (7) The instantaneous power; (8) Time-averaged 
power; (9) Time-scale over which non-thermal particles are injected; (10) Ratio of the observation time to 
particle injection time-scale; (11) Deprojected distance. 
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not deviate far from the equipartition value, -Begw Both of these conditions are satisfied by 



the parameter sets as shown in Table 3.5 where the ratios B/Begu and Nnth/Nx are given. 

For each source, the time delay Ati2 between the ejection of the two blobs that form the 
knots are nearly equal to the time difference between the ejection of the first blobs of Knot 
A and Knot B, t^^. This gives an overall single time-scale of activity for each source which 
ranges from 10^^ — 10^^ s and can reproduce the knot properties, as had been assumed in 
the development of the simple internal shock model. This result is important since if it had 
not been true, a more complex temporal behavior would have to be proposed, wherein the 
jet structure is due to variability of the source at two different time-scales, the first being the 
time difference between the ejection of two blobs that form a knot, and the second being the 
time difference between the activities that produced the two knots. 



The cross-over time, Ton,k, is determined here by equation (3.13). Since, f32K — I^ik ~ 
(ri-r2)/r ^ 0.3, Toat.x ~ 10^^(Asx/5x lO^^cm) s. During this time, i.e., the time when 
there is injection of particles into the system, the knot would travel a distance ~ cTqn ~ 
cVToN,K ~ 50 (r/5)(Ax/c/5 x lO^^cm) kpc. This is a significant fraction of the total 



observable distance traveled by the knot, from formation ^ etc ~ 100 kpc (see Table 3.5 1 
to the termination of the jet in the radio lobe, ^ 200 kpc. Hence, as shown in Table |3. 5} it is 



possible to fit the spectrum of all the knots with an observation time to^K, which is less than 
the crossing-over time Ton,k, implying that there is continuous injection of particles into 
the system. In this scenario, the synchrotron and IC spectra will have a break corresponding 
to a Lorentz factor 7c where the cooling time-scale equals the observation time §3.1[ 



The knots of 3C371 and Knot A of 1136-135 are unique in this sample, since their the X- 
ray flux lies below the extrapolation of the radio-to-optical spectrum to X-ray wavelengths. 
This allows for the interpretation that the X-ray flux is due to synchrotron emission [[7811981 . 
For Knot A of 1136-135 and Knot B of 3C371, this implies that the spectral break for the 



synchrotron emission occurs at the X-ray regime (Figure 3.5), which in turn indicates that 
these sources are relatively younger. Indeed, the ratio of the observation time to the injection 



time, to,K /ToN,K, for these sources are smallest (Table 3.5 ). On the other hand, for Knot A 
of 3C37 1 , the spectral break can occur at the optical band even if the X-ray flux is interpreted 



as being due to synchrotron emission (Figure 3.5 1 and hence this source need not be relatively 
young. However, this is only possible if to,K < Ton,k and there is continuous injection of 
particles. Otherwise, a sharp cutoff in the spectrum at the optical band would have occurred 
and the X-ray emission would not be due to synchrotron emission. 
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Figure |3.5| shows that the radio, optical and X-ray spectral indices for different knots may 
vary and highlights the need for more spectral measurements in all bands. A definite pre- 
diction of this model is that for most knots, the X-ray spectral index should be equal to the 
radio spectral index, indicating that the X-ray flux is due IC/CMB. Such spectral constrains 
would be particularly important since, although it has been demonstrated here that the inter- 
nal shock model can explain the broadband spectrum of these sources, there could be other 
models that may be physically and observationally more favorable. Analysis of the Very 
Large array (VLA) and Hubble Space Telescope (HST) images of 3C273 have shown that 
the optical and radio spectral indices are different, indicating the presence of an additional 
emission mechanism for the source [I109II . Also, the X-ray flux may be due to a second 
population of non-thermal electrons, rather than being the IC/CMB spectrum of the same 
distribution that produces the radio and optical emission UllOH . An argument in support of 
this is that since in the IC/CMB model the X-ray emission is due to electrons that are only 
a factor ten more energetic than those which produce the radio, a source in which the X-ray 
flux falls rapidly from the center should also exhibit a similar decrease in radio emission. 
However this feature is not observed (e.g. 3C273 and 1354-1-195). Moreover, the jet power 
required in the IC/CMB model can be very large, ^ 10'^^ ergs/s, which may be larger than 
the power inferred from the giant radio lobes (< 10^^ ergs/s). While the former argument 
may not strictly be applicable to the internal shock model (since each knot is a separate entity 
and the distance from the source is not a measure of the age of the source), the power re- 
quirement for some sources may indeed be very large, for e.g. Knot A of 1 136-135 requires 



Pave ~ 2 X 10^^ ergs/s (Table |3.5[). However, the energy requirement may be decreased if 



the magnetic field is sub-equipartition e.g. 3C273 (Table 3.5 1. Thus it is desirable to obtain 
direct observational signatures, such as spectral indices, to discriminate between models. 

A realistic description of the knots is more complicated than the simple model considered 
here. For example, the forward and reverse shocks that should form when the blobs collide, 
may provide different injection rates and at different locations within the Knot. However, the 
physics of these shock formations and the subsequent acceleration of particles is complicated 
and unclear, especially if they are mediated by magnetic fields. In the future, results from 
sophisticated numerical simulations could be compared with higher resolution data (which 
can resolve the internal structure of the knots) to prove (or disprove) the internal shock 
model. 
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3.3 A Two Zone Model for the knots of M87 

The models discussed in earlier sections can explain the broadband spectrum from the knots 
of many AGN jets. However they fail in case of M87, a nearby giant elliptical galaxy, at a 
distance =16 Mpc possessing an one-sided jet with projected distance ^ 2 kpc. The jet 
is bright in radio, optical and X-ray energies. The jet structure is very well studied in radio, 
infrared, optical and X-ray energies ll49l[Tl^[TT3l[TT4l[TT5l[T00ll . The flux and the spectral 
indices at X-ray energies indicate a possible continuation of synchrotron emission of the 
radio-to-optical spectrum with the change in the spectral index beyond optical energies [|116[ 
I117L Simple theoretical models, namely the continuous injection model ir76llll8[ril9[ll20ll 
( §3.1| ) and the one-time injection model [|105[ 1761 II 061 were unable to explain the observed 
X-ray flux and/or the spectral index. The X-ray flux predicted by the continuous injected 
model is more than the observed flux whereas the one-time injection model with pitch angle 
scattering under predicts the X-ray flux and the one without pitch angle scattering fails to 
predict the observed X-ray spectral index HIOOB . 

We propose a two zone model to explain the non-thermal emission from the knots of M87 
jet. In a two zone model, particles are accelerated in a region, namely acceleration region 
(AR) (probably around a shock front), and subsequently cool off through radiative processes 
in an associated cooling region (CR). 

3.3.1 The Model 

We consider the acceleration of a power-law distribution of particles (which may be a relic 
of a past acceleration process) at a shock front and cooling via synchrotron radiation in a 
homogeneous magnetic field. We treat the present scenario as two zones: one around the 
shock front where the particles are accelerated (AR) and the downstream region where they 
lose most of their energy through the synchrotron process (CR). This model is then used to 
explain the radio-optical-X-ray spectrum of the knots in the M87 jet. We assume the CR 
to be a spherical blob of radius R with tangled magnetic field Bcr and the AR is assumed 
to be a very thin region with magnetic field Bar- A power-law distribution of electrons is 
continuously injected into the AR characterized by an acceleration time-scale tacc- Particles 
are then accelerated at a rate I /tacc to a maximum energy determined by the loss processes. 
The AR is assumed to be compact and the emission from the CR mainly contributes the 
overall photon spectrum. 
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Fig. 3.6: Grey-scale image of M87 in radio (top), optical (middle) and X-ray (bottom). Figure repro- 
duced from Wilson & Yang llT2ll . 
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The kinetic equation governing the evolution of electrons in AR is given by 



c?n(7,t) d 
dt ^7 



CarI 



7 



n(7,t) 



n(7,t) 



+ Q(7) 



(3.22) 



where 



Q(7) d-f = qoj ^ d-f for -fmin < 1 < lb 



(3.23) 



Here 7 is the Lorentz factor of the electron, tesc is the escape time-scale and C,ar = c ^\r- 



Equation (3.22) can be solved analytically using Green's function and the electron distribu- 
tion for an energy-independent tacc and t^sc at time t is given by 111221 



n 



7 



7 



Imax 
-a 



a-1 



X 7„ 



dx 



(3.24) 



where a = tacc/tesc and the lower limit of integration Xo is given by 



Xq 



In 



+ 



7 lr> 



exp{t/tacc) 



(3.25) 



7« 



1 / {Car tacc) > 7b is the maximum Lorentz factor an electron can attain in AR. For 



t » tacc equation (3.25) can be approximated to be 7mm as the injection term in equation 



p.24| ) vanishes for x < •jmin- 

The evolution of the electrons in CR is governed by the equation 



dt 



87 



[CcRl''N{^,t)\+QAR{l) 



(3.26) 



where Ccr = c -^cr ^'^d the last term is the injection from AR, Qar{i) = n{'y)/te 
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Fort >> t„ 



7 



7n 



a-1 



MiN(^irib) 



X 





'1 1 








_ '~imax _ 



dx 



(3.27) 



The distribution of electron at time t in CR from equation (3.26 1 is given by 



iV(7,t) 



CcR 7^ 



QARix) dx 



(3.28) 



where To = 7/(1 - 7 Ccr i)- 



From equation (3.24 1, it can be shown that the injection into CR (for a + 1 > p) is a broken 
power-law with index — p for 7 < 7b and — (a + 1) for 7 > The synchrotron losses in 
CR introduces an additional break 7c in the electron spectrum depending upon the age of CR 



(W and the 5c;/? (^3.1). 



7c = — — (3.29) 

^CR f^obs 

The electron spectrum in CR at tobs can then have two different spectral shapes depending 
on the location of 7c with respect to 

(i) 7c > 7^: The final spectrum wiU have two breaks with indices 

7 ^) Imin ^ 7 ^ 7fe 

iV(7, Us) cc { 7-("+i), 76 < 7 < 7c (3.30) 

7c < 7 < 7max 
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(ii) 7c < 7fe: In this case the indices are 

7"^, 7mm < 7 < 7c 

iV(7,t„,,) oc <( 7-(p+i), 7c<7<76 (3-31) 

^—{a+2) ^ < 'y < 'y 
I ) /c ^ / ^ Jmax 

The resultant synchrotron emissivity can then be calculated by convolving A^(7, t) with 
single particle emissivity averaged over an isotropic distribution of pitch angles (equation 
( |247] )). 



The predicted spectrum depends on nine parameters, which are qo, a, jmin, 76, Jmax, P, 
BcR, R and tobs- Here a and p are estimated from the radio-to-optical and optical-to-X-ray 
spectral indices; qo and Bcr are constrained using the observed luminosity and equipartition 
magnetic field. For R we assume the physical sizes measured in radio [1123II . The age of the 
knot tobs is chosen to introduce a break in the observed spectrum at optical band and 7^ is 
fitted to reproduce the observed X-ray flux. 7^™ and 7max are used as free parameters and 
are fixed at 5 and 10®. 



3.3.2 Results and Discussion 



We applied the above model to explain the knots D, F, A and B of M87 jet (Figure 3.6 ). The 



results of the fitting are shown in Figure 3.7 and the parameters used for the fit are given in 
Table 3.6 The spectrum of knot E can be explained by simple continuous injection model 
and the parameters we quote corresponds to continuous injection model. We did not model 
knot C due to significant differences in X-ray-optical properties. The optical image of knot 
C is a diffuse region with a single maximum whereas in the X-ray image there exists two 
distinct maxima coincident with the diffuse optical knot HIOOL 



For all the fits shown in Figure 3.7 7c < 7f,. However, one can fit the spectrum with 7c > 
7b with proper choices of the parameters a, 75 and 7c. This degeneracy arises due to the 
unavailability of the ultraviolet (UV) spectral index since the present model predicts the 



corresponding particle spectral index as —{a + 1) or — (p + 1) (equations (3.30) and (3.31 1) 
depending upon the above two conditions. Future observations at these photon energies may 
help in validating the present model and also will remove this degeneracy. Also, to obtain 
a precise values for p and a, spectral indices at radio, optical and X-ray energies should 
be known accurately. It should be noted here that in general tacc and tgsc can be energy 
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Fig. 3.7: The spectral fit for the knots of M87 with the observed fluxes in radio, optical (Perlman et 
al. LI 24,1 and X-ray (Perlman & Wilson II 00.1 ) . Knot labels are the same as used in the 
literature ifTITTl . 
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Table 3.6: Model Parameters for the Knots of M87 



Knot 


9o 


a 


P 


76 


BcR 


tobs 


7c 


R 






(10-12) 




(10«) 


iW^ G) 


(10^ s) 


(105) 


(pc) 


D 


7.5 


1.75 


2.35 


1.8 


9.3 


1.6 


5.7 


12 


E 


4.7 




2.36 




5.9 


2.5 


9.0 


17 


F 


0.6 


2.3 


2.3 


2.1 


4.5 


4.7 


8.2 


29 


A 


1.0 


2.45 


2.29 


2.0 


4.7 


3.5 


10.1 


55 


B 


0.9 


2.75 


2.3 


1.4 


4.7 


3.9 


9.0 


50 



Columns: (1) Knot name; (2) Normalisation of the power-law injected into AR (for knot E it is the normali- 
sation of the power-law injected into CR (see text)); (3) Ratio between the acceleration time-scale and escape 
time-scale in AR; (4) Index of the power-law spectrum injected into AR; (5) Maximum energy of the electron 
Lorentz factor injected into AR; (6) CR magnetic field; (7) Time of observation; (8) Break Lorentz factor of 
the electrons due to synchrotron cooling in CR; (9) Size of CR measured in radio li,123.l . 
For all cases, the minimum Lorentz factor {jmin) injected into AR is 5 and maximum Lorentz factor {'^max) 
attained in AR is 10^. 



dependent. In such a situation the solution (equation (3.24)) may differ from its form and the 
index beyond 75 may not be the one discussed above. 

A possible scenario of the present model is where the AR is a region around an internal shock 
following an external shock. The electrons injected into the AR can be those that are already 
accelerated by an external shock and are advected downstream to be accelerated further by 
the internal shock [|125[|126l . Alternatively, re- acceleration of a power-law electron distribu- 
tion by turbulence at boundary shear layers can also be another possible scenario [1127111281 . 
Inclusion of these scenarios in their exact form into the present model will make it more 
complex and is beyond the scope of the present work. 

Perlman & Wilson (2005) proposed a modified continuous injection model where the vol- 
ume within which particle acceleration occurs is energy dependent HlOOl . This is expressed 
in terms of a filling factor face which is the ratio between the observed flux and the flux pre- 
dicted by the simple continuous injection model. They found face declining with increasing 
distance from the nucleus suggesting particle acceleration taking place in a larger fraction 
of the jet volume in the inner jet than the outer jet. The energy dependence of face also 
indicates that particle acceleration regions occupy a smaller fraction of the jet volume at 
higher energies. Even though the model is phenomenological, it indicates that the process 
of high-energy emission from the knots is as complicated as their physical region. However, 
the mechanism responsible for the filling factor is not explained. 

Liu & Shen (2007) proposed a two zone model to explain the observed spectrum of the 
knots of the M87 jet. In their model electrons are accelerated to relativistic energies in the 
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acceleration region (AR) and lose most of their energies in the cooling region (CR) through 
the synchrotron process. They considered that the AR and CR are spatially separated and 
introduced a break in the particle spectrum injected in the CR through the advection of 
particles from the AR to the CR. This along with the cooling break in the CR produces a 
double broken power-law with indices —p, —{p+ 1) and —{p + 2) which is then used to fit 
the observed spectra. However, the present model assumes that the AR and CR are cospatial, 
supporting a more physical scenario where electrons accelerated by the shock cool in its 
vicinity. 

We explored the possibility of the present model to reproduce the X-ray flux of other FR 
I galaxies (detected by Chandra) which are observed to have lower radio luminosity and 
relatively smaller jets when compared with FR II galaxies. The X-ray emission from the 
knots and/or the jets of the FR I galaxies, namely 3C 66B [fT29l . 3C 346 [fBOll . Cen A ifBTll 
and 3C 296 II132II , listed in the online catalogue of extragalactic X-ray jets XJET^, which are 
not explained by synchrotron emission from simple one zone models, can be reproduced by 
the present model. 



^http://hea-www.harvard.edu/XJET/ 



Chapter 4 



Boundary Shear Acceleration in the Jet 
ofMKN 501 



MKN 501 is a nearby BL Lac object (z = 0.034) and also the second extra galactic source 
detected in TeV photon energies by ground-based Cherenkov Telescopes II133II . It was later 
detected in MeV photon energies by the satellite -based experiment EGRET [I134L Radio 
images of MKN 501 show a jet emerging from a bright nucleus [11351 1136[ 11371 11011 . 
The high-resolution (milliarcsecond) radio images show a transverse jet structure with the 
edges being brighter than the central spine commonly referred as "limb-brightened" structure 
Ill35[|136lll01ll (Figure 4.1 ). This feature is usually explained by the "spine- sheath" model 
where the velocity at the jet spine is larger compared to the velocity at the boundary. Such 
a radial stratification of velocity across the jet arises when jet moves through the ambient 
medium and the viscosity involved will cause a shear at the boundary. Three-dimensional 
hydrodynamic simulations of relativistic jets [I138II and two-dimensional simulations of rel- 
ativistic magnetized jets [11391 also support the presence of jet velocity stratification due to 
its interaction with the ambient medium. The existence of velocity shear at the jet bound- 
ary was first suggested by Owen, Hardee & Cornwell [|49l[ to explain the morphology of 
M87 jet. Perlman et al. [[140[[ later confirmed it through the polarization studies of M87 
jet. Limb-brightening can occur in a misaligned jet with velocity stratification. For a proper 
combination of flow velocities, one will see a Doppler-boosted image of the boundary com- 
pared to the less boosted spine due to relativistic effects [[141[[142| . A possible consequence 
of the velocity shear is the alignment of the magnetic field at the boundary parallel to the flow 
velocity due to stretching of the frozen-in field lines of the plasma [[143[[ . The polarization 
angle observed at the jet boundary of MKN 501 is perpendicular to the jet axis indicating 
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Fig. 4.1: Grey-scale image of MKN 501 jet in radio at 15 GHz (top left), 22 GHz (top right) and 1.6 
GHz (bottom left and bottom right). Figures reproduced from Giroletti et al. jlOlll . 



a parallel magnetic field [|144[|137ll . However, it should be noted here that the polarization 
angle at the jet spine indicates a perpendicular magnetic field, and this along with the paral- 
lel magnetic field at the jet boundary can be an outcome of a dynamically dominant toroidal 
magnetic field structure [|144[ 11451 1146II . The radial velocity stratification of the jet can intro- 
duce Kelvin-Helmholtz instability, and the stability of jets against this instability was studied 
by various authors [M M [BOl [BIl . 

Giroletti et al. HlOlll have studied the limb-brightened structure of MKN 501 jet considering 
the differential Doppler boosting at the jet spine and the boundary [|142[|1411 . and concluded 
that the viewing angle (angle between the jet and the line of sight of the observer) of the 
radio jet should be more than 15°. However, high-energy studies of MKN 501 demand that 
the viewing angle of the jet should be ~ 5° in order to explain the observed rapid variability 
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and the high-energy emission [|152[ I153L Considering the fact that the 7-ray emission is 
originated from the inner part of the jet close to nucleus, Giroletti et al. HlOlll suggested that 
a bending of the jet may happen immediately after the 7-ray zone to explain the required 
large viewing angle of the radio jet. However, the mechanism required to bend the jets 
is still not well understood (jets deflected due to the pressure gradient in external medium 
are studied by Canto & Raga [43J; Raga & Canto Mendoza & Longair 1115411 )). and 
moreover the observed large bending of the jet in the radio maps can be the apparent one 
because of projection effects. This projection effects are even amplified when the jet is close 
to the line of sight. However it needs to be noted here that jets with large bending angle are 
indeed observed flUJ. 

The limb-brightened structure can also be explained if we consider the synchrotron emis- 
sion from the particles accelerated at the boundary, and this inference does not require large 
viewing angle. Eilek [11551 |48ll considered the acceleration of particles due to turbulence 
initiated by Kelvin-Helmholtz and Rayleigh-Taylor instabilities at the jet boundary. Parti- 
cles at the boundary can also be accelerated via shear acceleration [|46l l47l . and this case 
is considered in the present work. The acceleration of particles in a shear flow or by tur- 
bulence is well studied by various authors for both relativistic and non-relativistic cases 

[BTl [IM [Bll [IM [MD [El . 



4.1 Shear acceleration at MKN 501 jet boundary 



The particle acceleration process at the jet boundary can be described by the diffusion equa- 
tion in momentum space. The evolution of an isotropic phase space distribution is given by 

Eel 



dt 



p2 Qp 



dp J 



(4.1) 



where D(p) is the momentum diffusion coefficient. The characteristic acceleration time- 
scale can be written as 



P 



'd_ 
dp 



{p'D{p)) 



1 -1 



(4.2) 



If we consider a sheared flow, the electrons are scattered across different velocity layers by 
turbulent structures which are embedded in the shear flow. Berezhko [46] showed in such 
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case that there will be a net gain of energy in the electrons getting scattered, and this process 



is referred to as shear acceleration (^2.1.3). The momentum diffusion coefficient in case of 



a shear flow can be written as [|160[|164ll 

Ds{p) = xp'r (4.3) 

where r is the mean scattering time given by r ~ A/c with A being the mean free path and 
X is the shear coefficient given for a relativistic flow as HI 6511 

X = T^T^TTT-o TT ^ (4.4) 



15(r(r)2 - 1) \dr 

where r(r) is the bulk Lorentz factor of the flow and r is the radial coordinate of the jet cross 



section. Using equation (4.3 1, the shear acceleration time-scale (tacc,s) for t = TqP^ will be 

tacc,s = T~. I 7T (4.5) 

In case of turbulent acceleration(stochastic), the particles are scattered off by randomly mov- 
ing scattering centres and gets energized by a second-order Fermi acceleration. The momen- 
tum diffusion coefficient in this case can be approximated as HI 6411 

DM^f^{^)^ (4.6) 
where the Alfven velocityCV^) is given by 

Here B is the magnetic field and p the mass density of the jet. Hence, the turbulent acceler- 
ation time-scale(tacc,t) will be 

tacc,t = T^(^\ (4.8) 



(4-0 \Va/ 

For shear acceleration to be dominant over turbulent acceleration tacc,s < tacc,t- If we con- 
sider Bohm diffusion {(^ = 1) then the mean free path of the electron aligned to the magnetic 
field (A||) scales as the gyro radius (vg) [I166II . Ay ~ rj ^^'^'^ , where 77 is a numerical factor 
(?7 > 1 for magnetized particles) and 1) is the Lorentz factor of the scattered elec- 
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tron. Since the magnetic field at the jet boundary of MKN 501 is parallel to the jet axis (or 
toroidal) HI 37[ 11441 1145II . we consider r ~ Ay/c. Also if we consider 



or 

dr 



AT 

Ar 



(4.9) 



where AF is the difference between the bulk Lorentz factor at the jet spine and the jet bound- 
ary and Ar is the thickness of the shear layer, then the condition for shear acceleration to be 
dominant over turbulent acceleration will be 



Ar < 



47r p 



3(r(r) 



(4.10) 



If we consider the mass density of the jet is dominated by cold protons and if the number 
of protons is equal to the number of non-thermal electrons, then the jet mass density can be 
written in terms of equipartition magnetic field(_Beq) as 



mpB%{2a-l) 

167r irie C2 a ^rnin 



(4.11) 



and equation (4.10) will be 



Ar < 0.29 



r/7c2(Ar) 



eq 



nie nip (2 a — 1) 
a 7mm (r(r)2 - 1) 



(4.12) 



where a is the observed photon spectral index, nip is the proton mass and jmin is the 
Lorentz factor of electron responsible for the minimum observed photon frequency z/^m- 
The equipartition magnetic field can be expressed in terms of observed quantities as 



^eq — 



9.62 



1 



VaTi2a-l] 



G 



(4.13) 



where F{vmin) is the flux at the minimum observed frequency i^mm, di is the luminosity 
distance, V is the volume of the emission region and ax is Thomson cross section. Hence, 
for r(r)^ ^ 1 and a ~ 0.7, shear acceleration will dominate the particle spectrum at the jet 
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boundary of MKN 501 if the thickness of the shear layer 



Ar < 7.22 X 10^^ x 



^ I^F(lOMHz) 



910 mJy 







vlO/ 


[lo 


-5 




14 / 


R 



^obs \ 



10 J V1.6 GHz/ VlO MHz/ 



-V 



X 



1.5 pc 



15 
14 



pc 



(4.14) 



where R is the radius of the spherical region considered. (We assume 10 MHz as minimum 
observed frequency, and the flux at 10 MHz is obtained from the flux at 1.6 GHz considering 



the same spectral index. The flux at 1.6 GHz and R in equation (4.14) are obtained from 
a region around R.A 10 mas and declination —10 mas from Figure 4.1 (bottom left)). The 
corresponding equipartition magnetic field B(,q for T = 5 is 1.2 x 10^^ G. 

The electrons accelerated by shear acceleration cool via synchrotron radiation. The cooling 
time for synchrotron loss is given by 



cool 



Using equations ( 4.5| ) and ( |4.15[ ), we find 



'"""''^ ^ 1.5 X 10-12 



t. 



cool 



7 



B 



(4.15) 



X 



Ar 



10-9 pc 



1.2 X 10"3 G 

2 



H 

10 



Tir) 



X 



Ary 



(4.16) 



and since tacc,s ^ icoou shear acceleration dominates over synchrotron cooling. It can be 



noted that equation (4.16) is independent of the electron energy and hence the maximum 
energy of the electron will be decided by the loss processes other than synchrotron loss 
(which are not considered in this simplistic treatment). 

If we maintain the general form of mean scattering time r = tqp^, then for shear acceleration 
to dominate over turbulent acceleration the thickness of the shear layer (Ar) should be 



Ar < 1.7 X 10^ 



Top? (Ar) 



rfri 



(4 + 0(2a- 
_ a(4-07« 



cm 



It can be noted that equations (4.10) and (4.17) are equal, if we set in the latter ^ 

TqP^ = 7] Tg/C. 



(4.17) 



1 and 
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4.2 Particle Diffusion at the jet boundary and 
Limb-brightening 

Particles accelerated at the shear layer of the jet boundary diffuse into the jet medium before 
getting cooled off via synchrotron radiation. As the magnetic field at the jet boundary is 
parallel to the jet axis (or toroidal) [|137[ I144[ I145II . the radial diffusion of the electron into 
the jet medium is determined by cross-field diffusion. The cross-field diffusion coefficient 
can be approximated as [|167[ 11681 II 6611 



1 

3 r] 



TgC 



(4.18) 



where ri{> 1) is the scaling factor determining the field-aligned mean free path (see ^:4.1 1. 



The radial distance Rdiff that the electron diffuse before getting cooled can then be approx- 
imated as 



R. 



diff 



cool 



(4.19) 



Using equations (4.15) and (4.18) and considering the equipartition magnetic field, we get 



Rdi 



diff 



2.9 X 10" 



10 



B 



1.2 X 10-3 G 



pc 



(4.20) 



Since the thickness of the shear layer Ar <^ Rdiff (refer equations (4.14) and (4.20)), the 



thickness of the limb-brightened structure will be ^ Rdiff- This corresponds to an angular 
distance of 4.7 x 10""^ mas which is beyond the resolution of present-day telescopes. 

For r = tq p^, the cross-field diffusion coefficient will be 

1 



6Tq ^ 



(4.21) 



Using equations (4.15) and (4.19) we get 

Rdiff ^^.l^iW^ B-'tq-" p^ cm 



(4.22) 



and hence the thickness of the hmb brightened structure will be energy dependent for ^ 1. 
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4.3 Spectral index 



If we add mono-energetic particle injection term iS{p—po)) and particle escape term (— 1/tesc) 



in equation (4.1 1, then the steady state equation in case of shear acceleration for p > po and 



,^ = 1 can be written as 



dfs fs 



dp^ dp X To te 

and in case of turbulent acceleration it will be 

d'ft , . dft ft 



(4.23) 



dp"^ dp Tp te 







(4.24) 



1/2 

where ijj — — - 



g^a:^. If we substitute p = 1/x in equation (4.23 1 we get 

d^fs , dfs fs 



X 



_34JL _ 







(4.25) 



dx'^ dx "X, '^o ^e„„ 

Equations (4.24) and (4.25 ) can be solved analytically P 169I . The solutions are complex and 
are given by 



fs 



1 



X 



XToPte 



tts Ji I 2i 



X 



XToPte 



+ bsYj 2i 



XToPte 



(4.26) 



and 



P 



at J2 2i 



P 



+ htYA2i 



P 



(4.27) 



where Jn{z) and Yn{z) are the Bessel functions of first and second kind and a^, bs, at and 
bt are constants. For negligible escape (tgsc 00), using the limiting forms of Bessel func- 



tions ni70B . the solutions, equations (4.26) and (4.27), approach a power-law fs oc and 
ft oc p~^. The shear- accelerated particle number density will then be ns{p) oc p"^ and the 
corresponding synchrotron photon flux will be S^^shear oc z/^^/^. For turbulent acceleration, 
the number density will be independent of p {nt{p) oc p^) and hence the observed syn- 
chrotron photon flux will be a flat one S^,turb oc v^l"^ [I106L The spectral index map of MKN 
501 jet indicates a steep photon spectrum at the boundary and a flat spectrum at the spine 
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Fig. 4.2: Low resolution 1.6 GHz-4.8 GHz spectral index map of MKN 501 jet. Figure reproduced 
from Giroletti et al. ifToTI . 



mOllK Figure 4.2 1. Hence, it can be argued that the shear acceleration may be dominant at the 
jet boundary of MKN 501 and turbulent acceleration at the jet spine. However, ^ is usually 
related to the turbulent spectral index [I171II which may be different at the jet boundary and 
jet spine. 



4.4 Discussion 

As the AGN jet moves through the ambient medium, the viscosity involved will cause a shear 
at the jet boundary, and hence acceleration of particles in these shear layers is unavoidable. If 



the shear gradient dT/dr is very steep or if the shear layer is very thin (equation (4.14)), then 
shear acceleration can dominate over the turbulent acceleration initiated by the instabilities 
at the jet boundary ll48l . Turbulent acceleration may play an important role at the interior 
regions of the jet HI 621 and can provide an alternative to explain the emission from the inter 
knot regions of AGN jets [1172111731 . The observed hard spectrum at the jet spine HlOlll also 
support this inference since turbulent acceleration can produce a hard particle spectra HI 6211 



(also shown in ^4.3). The electrons accelerated by the turbulence can be reaccelerated by 



shocks and can form a broken power-law electron spectrum. This can possibly explain the 



break in the radio-to-X-ray spectrum of the knots of FRI jets (^: 3.3 1. 



Chapter 4. Boundary Shear Acceleration in the Jet of MKN 501 



87 



Giroletti et al. HlOU calculated the jet viewing angle (9) using the correlation between the 
core power and the total power [|174| . They estimated the jet viewing angle to be within 
10° < < 27° by comparing the observed core radio power and the expected intrinsic core 
power derived from the correlation. However, this estimation may vary if the core flux den- 
sity variability is more than a factor of 2. Also considering the variation of the parameter 
values in the correlation with increased number of samples, this may not provide a strong 
constrain on the jet viewing angle. The estimate of 9 based on the adiabatically expand- 
ing relativistic jet model [I175II may not be a strong constraint as it considers a simplified 
situation. Also, the constrain is less severe in case of perpendicular magnetic fields, and 
observed polarization studies have indicated the presence of perpendicular magnetic fields 
at jet spine Ill44[|137l l. Stawarz & Ostrowski II159II proposed a model similar to the present 
one; however, their aim was to show the observational implications of the two-component 
particle spectrum (power-law distribution with high-energy pile-up) formed at the boundary 
shear layer and the complex beaming pattern. 



Chapter 5 



A Two zone Model for Blazar Emission 
Mechanism 



The spectral energy distribution(SED) of blazars are characterized by a typical double-hump 
feature. The first component peaks at IR/optical energies for low energy peaked blazars 
and at UV/soft X-ray energies for high energy peaked blazars. Whereas, the second com- 
ponent peaks at hard X-ray/7-ray region. The first component is generally modelled as the 
synchrotron emission due to the cooling of relativistic electrons in a magnetic field while 
the second component is considered to be produced due to the inverse Compton scattering of 
soft photons by the relativistic electrons. The target photons for the inverse Compton process 
can be either synchrotron photons themselves (SSC) HI 76111771 II 02[[83l or external radiation 

(EC) (^233 [inHi 1691 [m. 

Besides their non-thermal continuum, blazars also exhibit rapid and strong variability (for 
example MKN 421 [[ISOlIISIlIIHl, MKN 501 [USSlIIll, PKS 2155-304 [[T85l[T86l[T87]l . 
etc.). In some cases the variability detected at different frequencies are non simultaneous 
and the flaring activity is referred as hard lag or soft lag depending upon the observed flare 
pattern. In a hard lag, the low frequency flare leads the high frequency one and for the soft lag 
it is the other way. Takahashi et al. H 17611 reported a soft lag in ASCA observation of MKN 
42 1 , while Fossati et al. UlSOi I181II reported a hard lag in the BeppoSAX observation of MKN 
421 in 1997-1998. Similarly for PKS 2155-304, a soft lag was reported by Chiappetti et al. 
II185II and Kataoka et al. II186II in BeppoSAX observation of 1997 and ASCA observation of 
1994 respectively. However, Takahashi et al. HI 8211 contradicted the conclusions of Fossati 
et al. [|180[ I181II by reporting no spectral lag for the same data set. Also, Edelson et al. 
II187II reported no spectral lag for PKS 2155-304 from XMM-Newton observation in 2000. 
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Nevertheless, the study of short time-scale variability at different wavelength region and 
their interrelations is very important to understand the geometrical and causal connections 
between emission regions and emission processes. 

The flux variability observed in blazars have been studied by several authors using one zone 
and two zone models HSOllHllIIOllHSl. Kirk et al. [[821 and Kusunose et al. [HI assumed a 
two zone model where particles are accelerated in a region presumably by a shock and escape 
into a cooling region where they lose their energy by radiative processes. On the contrary, 
one zone models assume instantaneous particle acceleration and do not consider the particle 
acceleration process separately. Chiaberge & Ghisellini HI 0211 explained the short time vari- 
ability observed in MKN 421 by dividing the emission region into thin slices. Moderski et 
al. [I188II modelled the variability detected in 3C279 considering SSC and EC processes. 

We consider a two zone model, namely acceleration region and cooling region, to study the 
spectral behaviour of blazars. This model is similar to the one discussed in §3.3| used to 
model the knots of M87 jet. However here we assume mono energetic particles are injected 
into the acceleration region which are then accelerated. We studied the effect of particle 
acceleration time-scale on the variability pattern as well as spectral evolution under two 
different scenarios: (a) acceleration time-scale is independent of particle energy and (b) 
acceleration time-scale depends on particle energy. 



5.1 Model 

We model the broadband emission of blazar from a spherical blob (cooling region) of radius 
R permeated by a tangled magnetic field B, moving down the jet with a bulk Lorentz factor 
r. The Doppler factor of the blob is given hy S = [T(l — (3 cos 6)Y^, where 6 is the angle 
made by the jet with the line of sight. 

Relativistic electrons are injected into the cooling region from an acceleration region, around 
a shock front moving with velocity Vs- The acceleration region is continuously fed with 
mono-energetic electrons with Lorentz factor 70 at a rate Qq electrons per unit volume. It is 
assumed that electrons are accelerated by diffusive shock acceleration process fTl\ and their 
evolution (5(7, t) is described by. 



dt d'y 



I acc 



d'-y 



syn 



Qii,t) 



Q{i,t) 



Qo'^(7-7o) (5.1) 
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where 



~dt 



7 



(5.2) 



and 



d'-y 
~dt 



4 aT 



syn 



3 mp c 



7 



R2 

Stt 



(5.3) 



are the particle acceleration rate and synchrotron loss rate of the particles in the acceleration 
region. Here tacc and tesc are the acceleration and escape time-scale of the electrons respec- 
tively and Bacc is the magnetic field in the acceleration region. In the frame work of diffusive 
shock acceleration, acceleration time-scale can be written as 



tn 



20 c 
3^ 



nif, c 
eB 



^7 



(5.4) 



where 



20 c 
3^ 



rrie c 
eB 



(5.5) 



The escape time-scale tgsc is assumed to be atacc where a is a constant factor. 



For energy independent scenario, 7 in equation ( 5.4 1 is replaced by a constant value, 7e/ / [183 



and for energy dependent scenario, equation (5.4) is maintained as such. The dependence of 
tesc on tacc is kept unchanged in both the cases. The analytical solutions of equation \5.l\ , 
under these two physical conditions, have already been discussed by ll82l . 



Accelerated particles are injected into the cooling region from the acceleration region at 
a rate Q{'j,t) /tesc where they their lose energy via synchrotron and SSC processes. The 
evolution of the particle spectrum A^(7, t) in the cooling region is governed by 



dNi^,t) 
dt 



d_ 



d'-y 
dt 



iV(7,t) 



loss 



+ 



N{i,t) g(7,t) 



T 



(5.6) 



where total energy loss rate of an electron. 



) loss V dt 



+ 



syn 



d'J 

~dt 
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The synchrotron loss rate is given by 



Itt 



syn 



4 a^^,Bl 
3 nifC Svr 



(5.8) 



and the energy loss rate due to SSC process is given by 



4 ClT 2rr 

7 Uph 



(5.9) 



Uph is the energy density of synchrotron photons in the cooling region and can be calculated 
from the synchrotron specific intensity. Here Bq is the magnetic field in cooling region and 
ax is the Thomson cross-section. 



The equation (5.6) is solved numerically using finite difference scheme HlOSi 11021 . The 
specific intensity of synchrotron radiation at frequency u is given by 



/.(^,t) = ^(l-e--«^) 



(5.10) 



where, synchrotron emissivity j„ and synchrotron self- absorption coefficient k^, are given by 



(equations (2.47) and (2.44)) 



"ymax 



^-W = 4^ / iV(7,t)P(7,z/)f/7 



(5.11) 



and 



^max 

1 /• iV(7,t) 



d 



TrmgZ/^ J ry (^rj2 _ dt 



[7(7'-l)^P(7,^)] 



(5.12) 



respectively. Here, P(7, u) is single particle emissivity (equation (2.36)). The synchrotron 
photon energy density is then calculated using 



Uph{7,t) 



Is{t^, t) du 



(5.13) 
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where lfT89]| 

i^maxil) = min[us^rnax,^rne / Ah'j] (5.14) 

Here Us^min and 

^ s.max arc the minimuiTi and maximum frequency of the synchrotron photons. 
The specific intensity of inverse Compton scattering in the Thomson regime is calculated by 
using the standard formulation [51 J. 

Different possible mechanisms for the generation of flare have been discussed in the liter- 
ature. For example, Kusunose et al. [83] simulated flares by changing relation between 
acceleration time-scale(tacc) and escape time-scale(tesc) in the acceleration region for certain 
duration. Mastichiadis & Kirk [|8TI discussed flares due to sudden changes in, (i) electron 
injection rate Qq\ (ii) maximum attainable energy of electrons '^max and (iii) cooling region 
magnetic field B. We simulated flares by increasing the injection into the acceleration region 
for a very short duration over and above the steady state injection. Finally, the observed flux 
is computed taking into account the Doppler boosting [|36ll and cosmological effects. 



5.2 Results and Discussion 
5.2.1 Spectral evolution 

We considered the quiescent spectrum of MKN 421 as a test case and deduce the best fit 
parameter set for both energy dependent tacc and energy independent t^cc- Values of the 



parameters are tabulated in Table 5.1 which are taken as the standard values and all the 
results presented in this work are based on these values. With these set of parameters tacc 
turned out to be 1.4 x 10^ s (in source frame) for energy independent case, while O.57 s 
(in source frame) for energy dependent case. Hence the particle acceleration rate ('j/tacc) 
becomes 0.70 x 10^^7 s^^ and 2 s^^ for energy independent and energy dependent cases 
respectively. 



The fitted spectra of MKN 421 for both the cases of tacc are shown in Figure 5.1 The 
slope of the spectrum before synchrotron peak is decided by a, which relates the escape 
and the acceleration time-scale of electrons in the acceleration region. The escape of the 
particles from the cooling region produces a break at z/ ^ 10^'' Hz. As the second hump 
of the spectrum is produced due to the self-synchrotron Compton process, the break in the 
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Table 5.1: Parameters used for fitting MKN 421 



Region 


Parameters 


Energy independent tacc 


Energy dependent tacc 


Acceleration region 


7o 


1 f\r\ 

1.00 


1 nn 

1.00 




Qo (10-^#cm-^ s-i) 


2.1 


2.1 




Bacc (G) 


0.137 


0.132 




Vs (in units of c) 


1 


0.17 




7e// 


10^ 








5 X 10^ 


5 X 10^ 




a 


1.70 


1.70 


Cooling region 


Bo (G) 


0.18 


0.18 




i?(cm) 


2.4 X IQis 


2.4 X 10^6 




r(s) 


2.5R 

c 


2.5R 

c 




S 


12 


12 



synchrotron component of the spectrum is also reflected there. 

As shown in Figure [5T[ the steady state spectrum of MKN 421 can be reproduced suc- 
cessfully by the two acceleration scenarios though rate of acceleration of particles in these 
scenarios are different. Hence the steady state spectrum of blazar cannot uniquely reflect the 
possible underlying electron acceleration mechanism and its dependencies on the electron 
energy. In case of energy independent tacc, the rate of electron acceleration depends on the 
electron energy and is much less than the electron acceleration rate for energy dependent tacc 
(which is a constant in this case) for a wide range of 7. Hence the evolution of the spectrum 
for energy independent tacc will be much slower compared to the energy dependent tacc case 
as shown in Figure This feature in the spectral evolution will also be reflected in the 
flare phenomena, if the flare is produced due to a sudden increase in the rate of injection of 
particles Qo in the acceleration region for a small duration. 

For energy independent case the system attains the steady state on the 8th day (in the ob- 
server's frame) of evolution. Although, in the energy dependent scenario, the spectrum at- 
tains the steady state faster but to make a comparative study of flares, we considered the 
steady state spectrum in both the scenarios only after the 8th day of evolution. We discuss 
below the characteristics of the flare phenomena under the two different scenarios of particle 
acceleration process discussed above. 
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10 12 14 16 18 20 22 24 26 28 

log V (Hz) 

Fig. 5.1: Two zone model fitting of MKN 421 archival data considering energy independent accel- 
eration time-scale (top) and energy dependent acceleration time-scale (bottom) using the 
parameters given in Table 15. 11 The observed data are taken from Kino et al. II1901 . 
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Fig. 5.2: Evolution of the model spectra in case of energy independent acceleration time-scale (top) 
and energy dependent acceleration time-scale (bottom). The spectra are plotted for every 
0.1 days up to 1.5 days and there after at 2nd, 3rd, 6th and 8th days (in observer's frame). 
The steady sate is arrived after 8 days in both the cases. 
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5.2.2 Flare 

As mentioned above, we simulate a flare by increasing the electron injection rate into the 
acceleration region by a factor of 50 for a duration of 0.2 days (in the observer's frame) after 
the system reaches the steady state . 

Energy independent tacc 



The simulated light curves for energy-independent tacc are shown in Figure 5.3 It is evident 



from top panel of Figure 5.3 that the high energy flares lag the low energy ones in time (~ 
kilo s) and it is also seen that the light curves are more asymmetric at lower energies. The lag 
in the light curves arise since the rate of acceleration of electrons in the acceleration region is 
proportional to energy of the electrons and the high energy electrons which are responsible 
for the high energy synchrotron emission, take longer time to attain the required energy. The 
asymmetry in the low energy light curves can be attributed to the longer synchrotron cooling 
time- scales of low energy electrons compared to the high energy ones. It can also be seen 
that the variability amplitude increases with frequency of emission. These features are in 
qualitative agreement with the BeppoSAX observations of MKN 421 [|180[|1811 . The lower 



panel of Figure 5.3 describes the simulated flare patterns for SSC component of emission. 
Since, the synchrotron photons are Compton boosted by the same population of electrons, 
the variability features in SSC component are qualitatively same as in the synchrotron com- 
ponent. 

The variation of the spectral index with respect to the flux (hysteresis loops) have also been 



studied at different energies and are shown in Figure 5.4 The sense of the loops is anti- 
clockwise representing the hard lag in the system. 

Energy dependent tacc 

In this case the acceleration time-scale is proportional to the energy of the electrons and 
hence the rate of acceleration is independent of electron energies (l/ta,o)- For the param- 



eters given in table 5.1 we find t^ o = 0.5 s (in the source frame). As the rate of particle 
acceleration is high and same for all electron energies, the flares are initiated almost im- 
mediately after the extra injection into the acceleration region occurs. This fact is evident 



from near-simultaneous flare patterns shown in Figure 5.5 The light curves at optical/UV 
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0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 



Time in observer's frame (days) 

Fig. 5.3: Light curves at different frequencies for energy independent acceleration time-scale during 
a simulated flare created by the enhanced injection into the acceleration region (see text). 
The vertical dashed line represents the duration of the enhanced particle injection in the 
acceleration region. The top panel corresponds to the emission due to synchrotron process 
and bottom panel corresponds to the SSC process. 
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Fig. 5.4: The hysteresis curves showing the variation of spectral index with respect to flux during a 
simulated flare for energy independent acceleration time-scale. The top panel corresponds to 
the emission due to synchrotron process and bottom panel corresponds to the SSC process. 
The sense of rotation of the loops is anti clockwise. 
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frequencies in the synchrotron component show prominent breaks in their rising part. Since, 
the rate of acceleration is high and same for all energies, electrons in acceleration region 
attain higher energy in shorter time-scale compared to the escape time-scale, (tesc)- These 
high energy particles enter cooling region and get cooled giving rise to high energy emission. 
Then they join the freshly injected lower energy electrons and contribute to the low energy 
emission as well. The difference in the rates of these contributions give rise to breaks in the 
light curves at low energy emissions. This phenomena gives rise to a tendency that the light 
curves at optical and UV bands peak later than the light curves in the soft and hard X-ray 
bands. Although, this is not very prominent in the simulated light curves, it is reflected in the 



clockwise sense of the hysteresis loop shown in Figure 5.6 The absence of such breaks at 
very low energy emissions (namely far-IR band and below) can be justified by the fact that 
the escape time-scale from the cooling region is smaller than the time required for the high 
energy electrons to cool to these energies. 

Flares follow the same features of asymmetry as in the case of energy independent tacc- 
Similar features are reflected in the flares of SSC component also. 

Such variations in the photon spectral index during blazar flare (clockwise and anti-clockwise) 
are observed for many sources [|176[ II 851 II 80111911 . 

Table 5.2: Values of r] for different frequencies 



Frequency 


V 


(Hz) 




10" 


2.05 ±0.02 


10" 


1.71 ±0.09 


10" 


1.99 ±0.16 


1024 


1.67 ±0.04 



Dependence of lag time-scale on tacc 

As shown in the previous section, the energy independent acceleration time- scale gives rise 
to a hard lag in the flare patterns for different radiation frequencies. The variation of lag 
time- scale for different shock velocities Vg is shown in Figure [577| The lag time- scales are 
calculated with respect to the flare at 10^^ Hz. It is found that the lag time-scale varies with 
Vg following a power-law (oc Vs~'^)- The fitted value of 7] for different frequencies are given 



in Table |5.2[ As tacc depends on the shock velocity Vs {tacc ~ Vs ^), so the lag time-scale is 
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Fig. 5.5: Light curves at different frequencies for energy dependent acceleration time-scale during 
a simulated flare created by the enhanced injection into the acceleration region (see text). 
The vertical dashed line represents the duration of the enhanced particle injection in the 
acceleration region. The top panel corresponds to the emission due to synchrotron process 
and bottom panel corresponds to the SSC process. 
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Fig. 5.6: The hysteresis curves showing the variation of spectral index with respect to flux during a 
simulated flare for energy dependent acceleration time-scale. The top panel corresponds to 
the emission due to synchrotron process and bottom panel corresponds to the SSC process. 
The sense of rotation of the loops is clockwise. 



Chapter 5. A Two zone Model for Blazar Emission Mechanism 



102 




0.6 0.7 0.8 0.9 

shock velocity Vg (in units of c) 

Fig. 5.7: Variation of lag time-scale with shock velocity in case of energy independent acceleration. 

Filled circles corresponds to 10^^ Hz, open circles corresponds to 10^^ Hz, filled triangles 
corresponds to 10^^ Hz and open triangles corresponds to 10^^ Hz. 
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Chapter 6 



Summary and Conclusions 



The complete and comprehensive understanding of the physics of AGN jet requires the 
study of different aspects of jet in different wavelength ranges. Such studies may give rise 
to a global picture of jets describing the dynamics and radiation emission processes self- 
consistently. Although this is a mammoth task and requires a long time, an attempt has been 
made in this thesis to study certain aspects of jet physics phenomenologically and to un- 
derstand the possible inter-connections between them. Here, we studied the dynamics and 
emission processes from knots in jets observed in radio-optical-X-ray bands. Similarly, we 
considered blazars where we studied the limb brightening effect in Mkn 501 observed in 
radio wave band. This feature has intimate connection with physics of jet flow and the parti- 
cle acceleration process at jet boundary. Our work has been further extended to include the 
study of flares in blazars in infra-red-to-TeV energy band. 

To model the emission from the knots we used the archival radio, optical and X-ray data. A 
continuous injection plasma model, where non-thermal relativistic electrons are injected into 
an expanding spherical region with a tangled magnetic field, is used to study the broadband 

emission from the knots. Injected relativistic electrons lose energy by synchrotron process in 
the tangled magnetic field of the knot and by up-scattering the cosmic microwave background 
photons. The expansion also introduces an adiabatic loss. For a given observation time, the 
particle distribution in the emission region will be a broken power-law with a break at an 
energy where the cooling time-scale is equal to the observation time. The resultant spectrum 
is fitted to the observed data to yield the parameters of the model. The parameters obtained 
from the spectral fitting are physically reasonable and they are used to obtain the kinetic 
powers of the jets. 
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Above work has been further extended to include the dynamics of knots. Assuming that 
the knots are produced due to the collision of matter shells ejected randomly in time from 
the central engine, the complete kinematics of the shells are used to study their location of 
collision and the energetics. Non-thermal relativistic electrons are produced in the shock 
generated due the collision of shells. These electrons emit radiation by synchrotron and 
inverse Compton process as described above. Apart from fitting the observed spectrum, the 
other main conclusion of this work is that the location of the knots can be reproduced from 
the physically acceptable choice of parameters in the kinematics of shell collision. Therefore 
internal shocks can be considered as one of the viable mechanism of knot generation. It is 
to be noted that the timescales obtained from the kinematics of shells and internal shock, 
convincingly support the continuous injection scenario. 

A two zone model was proposed to explain the emission from the knots of the M87 jet since 
simple models involving continuous injection/one-time injection of non-thermal particles 
failed to reproduce the observed X-ray flux and the spectral index. In the proposed model, 
we consider the injection of a power-law distribution of particles into an acceleration region 
where they are accelerated further. The particles then escape from the acceleration region 
into a cooling region where it lose energy mostly via synchrotron radiation. The particle 
distribution in the cooling region will be a double broken power-law with one break at energy 
corresponding to the cutoff energy of the initial injected power-law into the acceleration 
region and the next break at energy for which the cooling time- scale equals to the age of the 
knot. The observed radio-optical-X-ray spectrum from the knots in M87 jet are reproduced 
by the resultant synchrotron emission. In its simplest form, the model does not consider any 
specific acceleration process but assumes an energy independent acceleration time-scale. 
The model can successfully reproduce the broadband spectrum from the knots/jets of other 
FRI galaxies, namely 3C 66B, 3C 346 and 3C 296, which are not explained by synchrotron 
emission from simple one zone models. 

High-resolution radio maps of the jet of the BL Lac object MKN 501 shows a limb-brightened 
feature and an explanation of this feature based on the differential Doppler boosting of a strat- 
ified jet, requires large viewing angle (> 15°). The viewing angle constraints inferred from 
the high-energy 7-ray studies of the source is very small (~ 5°). Since the 7-ray emission 
originates from the inner region jet, close to the central engine, this model requires the jet to 
be bent to accommodate the viewing angle conflict. However, the observed limb-brightened 
structure of the MKN 501 jet can be explained if we consider the shear acceleration of par- 
ticles at the boundary due to velocity stratification and their diffusion into the jet medium. 
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This inference does not require a large viewing angle as demanded by the explanation based 
on differential Doppler boosting of the jet spine and boundary. We have shown that shear 
acceleration dominates over turbulent acceleration at the boundary if we consider thin shear 
layer or a sharp velocity gradient. Also for the estimated set of parameters, shear acceler- 
ation time-scale is much smaller than synchrotron cooling time-scale allowing acceleration 
of electrons to be possible. The thickness of the limb-brightened structure will be decided 
by the distance electrons have diffused into the jet medium before loosing its energy via 
synchrotron radiation. However the estimated thickness is beyond the resolution of present 
day telescopes. Simple analytical solution of the steady state diffusion equation consider- 
ing mono-energetic injection and particle escape indicates a steep particle spectra for the 
electrons accelerated at the shear layer in comparison with turbulent acceleration. The radio 
spectral index map of MKN 501 jet is also observed to have steep spectrum at the boundary 
supporting the presence of shear acceleration. 

The temporal behaviour of the blazar emission is studied under the framework of a two zone 
model. The spectral evolution has been examined for two different physical conditions of 
diffusive shock acceleration mechanism, namely energy independent acceleration time-scale 
and energy dependent acceleration time-scale. The model is applied on the BL Lac object, 
MKN 421, to study the implications on the flare characteristics for the above mentioned 
conditions. We found that in case of energy independent particle acceleration, the photon 
spectrum evolves at a slower rate compared to the energy dependent case though their steady 
state spectra are not differentiable. The flare patterns at different frequencies show a hard 
lag in the energy independent acceleration scenario while they are near simultaneous in the 
energy dependent scenario. Hence, the presence/absence of time lags in the flare pattern 
has direct bearing on the underlying particle acceleration mechanism in a blazar jet. Also 
the presence of a break in the rising part of high energy light curves in the case of near- 
simultaneous flares suggests that the acceleration time-scale may depend on particle energy. 
In the case of energy independent tacc, it is also shown that the time-lag between two given 
frequencies has a power-law dependence on the shock velocity. Hence, a simultaneous multi 
wavelength study of blazar variability with good time resolution may be useful to constrain 
the physical parameters of the blazar jet and may also reveal the nature of underlying particle 
acceleration process which is crucial in understanding the dynamics of the jet. 

The models described here to study the different aspects of jets from different AGN can be 
improved further to determine and constrain the model parameters unambiguously if better 
quality data are available. This requires truly simultaneous long-term observation of AGN 
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using both ground-based and space-based telescopes in different wavelength bands. Partic- 
ularly, radio observations with high spatial resolution are necessary to study the structure 
of knots as well as the boundary layers of the jets. This will shed further light on the jet 
dynamics, instabilities and the possible particle acceleration sites along the jet. 

To study the flux variability of blazars it is important to have simultaneous measurements 
of flaring events at different wavelengths. The high sensitivity timing study is particu- 
larly important at GeV-TeV energies where blazars show the fastest flux variations. More 
observations of blazars with presently operating ground-based Cherenkov telescopes like 
MAGIC, HESS, VERITAS telescopes are extremely important in this respect. The upcom- 
ing Cherenkov telescopes with higher sensitivity like MACE and CTA will of course enhance 
the quality of data. This will in turn help us to improve the theoretical models to have better 
understanding of AGN jets. 

In fact the present work of blazar can further be extended to study the effects of extragalactic 
background light (EBL) on the blazar spectrum and possible estimation of EBL. This is an 
important issue in blazar research at high energies and this will be pursued in future. 
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